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ABSTRACT 


The methodlogy of developing finite life distributional 
Goodman diagrams and surfaces is presented in this paper. The 
Goodman surface and diagram presents allowable combinations of 
alternating stress and mean stress to the design engineer. 

The combined stress condition presented in these surfaces and 
diagrams is that of an alternating bending stress and a constant 
shear stress. The finite life Goodman diagrams and surfaces 
are created from strength distributions developed at various 
ratios of alternating to mean stress at particular cycle life 
values. 

The conclusions drawn in this report indicate that the 

4 

Von-Mises Hencky ellipse , for cycle life values above 10 
cycles, is an adequate model of the finite life Goodman dia- 
gram. In addition, suggestions are made which reduce the 
number of experimental data points required in a fatigue data 
acquisition program. 
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CHAPTER I 


INTRODUCTION 

Presently in the United States and abroad fatigue data, 
to a large extent, is presented in the form of the conventional 
S-N diagram. The S-N diagram's purpose is to graphically pre- 
sent strength data as a function of cycle life. This is done 
by testing a few specimens to failure at incremented stress 
levels. The data is plotted on a graph of log stress versus 
log cycles , (see Figure 1.1). This method of presenting fatigue 
data does not take into account one of the fundamental and most 
important fatigue aspects: the variability of the fatigue 

mechanism; i.e., even high quality test specimens, subjected 
to tightly controlled test conditions will rarely, if ever, 
fail at precisely the same cycle of life at a given stress 
level . 

Recently the American Society for Testing and Materials 
(ASTM) has suggested that a statistically significant number 
of specimens at each stress level be tested in order that a 
failure distribution be developed for each stress level (1, 
p. 9). The cycles-to-failure distributions which are developed 
can be used to construct a statistical S-N diagram. This 
particular diagram, as pictured in Figure 1.2, has a mean line 
as well as plus and minus three sigma lines. 


1 



Fig. 1.1 Conventional Alternating Stress Cycles To Failure Diagram. 
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It is possible through such statistical techniques as 
the use of probability plotting paper, the Chi-Squared and 
Kolnogorov-Smirnoff goodness of fit tests, and the computation 
of the four statistical moments to determine the best failure 
distribution probability density function. The probability 
density functions usually considered in such an analysis 
are the normal, lognormal, and Weibull distributions. The 
Weibull, because it is a three parameter distribution can take 
on many different shapes, from the exponential to the lognormal, 
by varying the three parameters. The Weibull probability den- 
sity function, because of this flexibility is the most flexible 
of the three distributions mentioned. 

The concept of failure distributions and strength distri- 
butions should be discussed briefly in order to avoid confusion 
in the later sections of this report. A failure distribution 
is derived directly from cycle to failure data. At a given 
stress level specimens will fail at particular values of cycle 
life. Even under the tightest controlled test conditions there 
will be variability in the cycle life of the individual test 
specimens. The failure distributions represent this variability 
in the test specimen's cycle to failure data. It is possible 
to form a histogram of the test specimen's cycle to failure data 
from which a failure distribution can be determined which ade- 
quately described the data. An example of such a histogram 
and distribution is given in Figure 1.3. Specifying the type 


Cycles to Failure 


Fig. 1 


3 An example of a Cycle to Failure Histogram and Distribution. 


* . 


6 


of distribution, normal, lognormal or Weibull, and the distri- 
bute n parameters, uniquely describes the cycles-to-failure 
data. 

A strength distribution is similar to a failure distribu- 
tion in that it is described by a distributional type and the 
corresponding distributional parameters. When a specimen fails 
it means that the stress has exceeded the value of the strength 
of the specimen. Cycles-to-f allure data reveals the percentage 
of specimens which have a strength less than the applied stress. 
It is possible to transform the cycles-to-failure distributions 
to strength distributions, A strength distribution describes 
the variability of the strength of a specimen at a specific 
value of cycle life and may be derived from a cycles to failure 
distribution. 

The discussion of transforming cycles-to-failure distribu- 
tions to strength distributions should be preceded by a discus- 
sion of some of the basic assumptions and restrictions which 
govern the generation of meaningful failure distribution dava. 

Of primary concern is that all test specimens be uniform in 
geometry and metallurgical properties. This is important since 
subsequent calculations of the strength distributions will as- 
sume that specimens tested at various stress levels came from 
a homogeneous population. In addition the number of test 
specimens should be dependent upon the variability of the data 
generated. The statistical significance of the desired data 
is dependent upon the number of test specimens which are run. 
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the more data points generated the nore positive the experi- 
menter can be of his theoretical distribution. It has been 
found that at lower stress levels the variability of the 
data increases which indicates that an increase in test 
specimens is warranted at these lower stress levels. 

The strength distribution may be oriented along the stress 
ratio axis of a dii ^ributional Goodman diagram. The ordinate 
axis of such a diagram presents values of alternating stress 
while the abscissa presents mean stress values. The strength 
distributions which are placed along axes where the ratio of 
mean stress to alternating stress is a constant form distribu- 
tional surfaces. These surfaces are formed when the strength 
distributions are connected by a mean line and plus and minus 
three sigma limit lines, (see Figure 1,4). The Goodman surfaces 
are of vast importance in reliability engineering where the 
interference of a stress distribution with the corresponding 
strength distribution is used to calculate the designed-in 
reliability of a mechanical part. 

The objective of this report is to clearly present the 
methodology of generating distributional Goodman diagrams. 

The accomplishment of this objective requires that the follow- 
ing subject areas be investigated: Chapters II and III explain 

methods of converting cycles to failure data to strength dis- 
tributions; Chapters IV and V develop methods of generating 
finite life Goodman diagrams and surfaces; Chapter VII directs 
it's attention to resolving the question of static strength 


Alternating Stress 



Fi*g. 1.4 Goodman Surface Formed By Strength Distributions 
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distributions to be used on the finite life Goodman diagram. In 
Chapter VIII empirical mathematical models of the Goodman diagram 
are discussed. Chapter IX explains the theoretical strength 
theories associated with the combined stress condition of alter- 
nating bending stress and mean shear stress. Chapter X recom- 
mends an empirical mathematical model of the Goodman diagram and 
an associated theoretical strength theory. Chapter XII suggests 
two methods of reducing the amount of experimental data needed 
to generate finite life Goodman diagrams as well as methods of 
obtaining cycles to failure distributions from these Goodman 
diagrams with a minimum amount of actual fatigue testing. 

The discussion of these subject areas requires that actual 
fatigue data be used in support of this effort. This investigator 
was extremely fortunate to have access to the complex fatigue 
data generated under National Aeronautical and Space Administration 
Grant No. 03-002-044 at The University of Arizona under the 
direction of Dr. Dimitri B. Kececiglu. The Combined stress con- 
dition under which this data was generated was that of an alternating 
bending stress and a constant, mean shear stress. Although the 
discussions in this report often apply themselves to this data, the 
concepts presented are applicable to the area of combined bending 
and shear stresses and in general to the broader area of any com- 
bined stresses in fatigue. 


CHAPTER II 


NASA METHOD OF GENERATING POLYGONS MW DISTRIBUTIONS 
2.1 Theory 

Mr. Richard E. Smith, Aerospace and Mechanical Engineering 
Department, The University of Arizona, in his master's report 
of August 1965, used data from Dr. H. T. Corten, Department 
of Theoretical ond Applied Mechanics, University of Illinois, 
to present a method of transforming theoretical cycle to 
failure distributions to cumulative strength polygons. The 
following is a summary of the methodology of that effort 
(2, p. III).. 

The theoretical cycle to failure distribution is first 
determined. Because goodness of fit tests allow only for the 
rejection of a distribution, it is possible that more than one 
of the three major fatigue probability density functions, normal, 
lognormal, and Weibull, will be accepted. It then becomes nec- 
essary for the investigator to choose which one fits the data 
best, at all of the various stress levels. Once this has been 
determined, the failure data is uniquely described by the 
theoretical failure distribution probability density function 
rather than the failure histogram of the sample data. This 
probability density function is symbolized by f(x). At each 
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of the j stress levels pictured in Figure 2.1.1, there is a 
theoretical normal, lognormal, or Weibull distribution which 
represents the failure data. 

The cumulative failure probability/ up to the x’th cycle, 
for the j stress level is given by 

x 

F (x) = Jf(x) dx (2.1.1) 

3 0 

Of interest to the designer is the strength distribution for a 
specified cycle life. Once this cycle life is specified a 
series of cumulative failure distributions can be calculated 
for each stress level. In essence, the calculation of F.(x) 
is equal to calculating the area bound by the theoretical 
failure probability density function and the cycle lino N 
(see Figure 2.1.2). Note the cumulative area in percent for 
each failure distribution to the right of the graph. Here it 
is important that the failure distributions are known for the 
full strength range in order that a cumulative failure probabi- 
lity of zero to one hundred percent is obtained (2, p. 23^. 

The physical significance of the cumulative failure distri- 
bution is the percent of specimens in which the stress has ex- 
ceeded the strength. The percentage of specimens at a given 
stress and cycle life with a strength equal to or less than the 
stress is new a known quantity. A plot, as shown in Figure 
2.1.3 can then be made of the cumulative failures in percent 
versus stress level, and is known as the cumulative strength 





Cumulative Failures *•“ Stress Cp. 
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histogram. From this emulative strength histogram it is possible 
to calculate the strength frequency histogram. The strength 
frequency histogram for the i'th cycle life is given by 

V s) = (p j+i (x) “ F j (x)) 

where F^. and F_. + ^ are given by Equation 2.1.1 and where there 
are H different stress levels to be considered. Since the 
above calculation for the strength frequency histogram is based 
upon inferences drawn from several stress levels it is impera- 
tive that the specimens used in all stress levels be from the 
sane statistical population and that uniform test conditions 
are maintained from stress level to stress level. 

Upon the determination of the strength frequency histogram, 
f.(s), the statistical operation for' goodness of fit, in this 
case the Chi-Squared test, can be conducted to insure that a 
theoretical distribution can be fitted to the histogram. For 
specific cycle of life values a theoretical strength distribu- 
tion is specified by one of the three theoretical distributions, 
specifically either the normal, lognormal or Ueibull distribu- 
tion. It is of major importance to note that this method is 
capable of determining the strength distribution in either the fatigue 
life or infinite life portion of the S-N curve (refer to Figure 2.1.1) 

This makes it unnecessary to conduct a Probit analysis 
or staircase test to determine the strength distribution (2, 
p. 31). 
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Beyond the endurance limit, area to the right of the knee 
where the S-N curve is horizontal, the failure distributions are 
independent of time, previously expressed in cycles. However, 
ot some stress levels it will be impossible to achieve a cumu- 
lative failures of one hundred percent because the test will 
be terminated at a pre -determined time before the specimen has 
failed. Fortunately the cumulative failure .distribution is 
known up until the test termination which allows one to calcu- 
late the strength distribution in an identical manner as 
previously discussed (2, p. 31). 

The result of the strength distribution calculations allows 
the construction of a statistical S-N diagram, as illustrated 
in Figure 2.1.4. 

An adequate range of stress levels is necessary so that 
a complete failure histogram from zero to one hundred percent 
can be developed. Within this range a sufficient number of 
failure distributions must be known at different stress levels 
so that there will be enough class intervals, governed by 
Sturges 1 rule, in the strength histogram. Time and economic 
considerations limit the number of stress levels, and conse- 
quently the number of failure distributions generated to from 
five to eight such levels and distributions. The use of this 
limited number of distributions would not give enough class 
intervals for an accurate strength distribution calculation. 

It is necessary to develop a digital computer program which 
will interpolate many stress-to-failure distributions from the 
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United number available from experimental data (2, p. 27). 

It is possible to have this program accomplish these interpola- 
tions and perform the methodology discussed for the transforma*' . 
tion of cycles to failure distributions to strength distribu- 
tions as well as accomplish all required statistical operations. 

2.2 Computer Method 

IWo computer programs were developed in order that the 
strength distributions could be calculated. The first reduces 
the cycles to failure data to failure distributions at each 
stress level. The second computer program transforms the failure 
distributions to strength distributions. 

Because of the large amount of fatigue data generated/ a 
digital computer program, in Fortran language, was used to 
reduce the failure data at various stress levels to failure 
distributions. Failure data in terms of stress level and cycle 
life is read into the computer. The computer is capable of 
calculating the following parameters: mean, standard deviation, 

coefficients of kurtosis and skewness, as well as of performing 
a Chi -Squared goodness of fit test. The computer program will 
use the normal distribution approximation in calculating the 
expected frequency in the Chi-Squared test when the sample 
data points are greater than thii'ty and will use the Student-t 
distribution when the sample number of data points is less than 
thirty. A flow chart, variable definitions and computer listing of 
the program are given in Appendix A. 
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The second computer program is used to determine the 
strength distributions from the failure distributions and is 
also written in Fortran language . The main steps of the pro- 
gram are listed below. A flow chart, cariable definitions and 
computer list of the program are given in Appendix B. 

1. For each of the experimental stress levels read 
into the computer the actual cycles to failure 
distribution parameters are developed. 

2. The computer will then calculate failure distribution 
parameters for interpolated stress levels at 200 

increments by straight line interpolation. 

3. The computer will then calculate the cumulative 
failure distributions for each stress level 
and a given series of log cycle life values. 

4. Next the computer will calculate the strength 
frequency histogram for each cycle life. 

5. The computer program then calls the computer to cal- 
culate theoretical distribution parameters from the 
strength frequency histograms based on a normal or 
lognormal distribution. These parameters include 
mean, standard deviation, coefficients of skewness 
and kurtosis. A goodness of fit tost, using the Chi- 
Squared test will then be performed to determine 
which of the normal of lognormal distributions fits 


the data best. 
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6. The parameters mentioned in five are then printed out in 
addition to the Chi-Squared values for each cycle life. 

This print out allows the investigators to determine 
which of the two distributions is a better fit. 

2.3 Results 

A complete description of Dr. Corten's testing progarm is found 
in Chapter V of Richard Smith's report, together with the analysis by 
Smith (2). Results of his analysis are summarized herein. 

The Chi-Squared goodness of fit tests indicated that for aluminum 
specimens the cycles to failure distributions more closely fit a lognormal 
distribution; whereas distributions of steel specimens fit either the 
notmal and lognormal distributions equally well (2, p. 76). It was 
observed that as the sample size was increased, the lognormal distribution 
fit the data better than the normal . 

The transformation of cycles to failure distributions to strength 
distributions was accomplished by assuming the failure data to be dis- 
tributed lognormal ly. The computer program previously described was 
used to determine the strength polygons for various cycles of life. The 
program output included the cumulative strength polygons, the mean, 
standard deviation, coefficient of skewness and kurtosis, and the Chi- 
Squared goodness of fit values for the strength distributions. 
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Analysis of this data indicates that the normal distribu- 
tion fit the strength data better than the lognormal distribu- 
tion. Coefficient of skewness is generally negative indicating 
a normal distribution. Coefficient of kurtosis value s, which 
should be 3.0 for the normal distribution fluctuate about a 
value of 3.0. In addition the Chi-Squared values indicated 
that for both type specimens the normal distribution represented 

V 

the data better than the lognormal distribution (2, p. 85). 

2.4 Discussion as to Validity 

The transformation of cycles-to-failure distributions to 
strength data which was proposed by Richard Smith.. was found 
to be appropriate and based on well founded principles. John 
Smith's methodology is identical in transforming cycles-to- 
failure data to strength distributions. Richard Smith applied 
the technique to Corten's data ( 3, p. 1X3) and John Smith to 


the KASA data. 


CHAPTER III 


GENERATION OF STRENGTH DISTRIBUTIONS 


3.1 Theory 

Smith (3) did extensive work in developing cycles to failure 
and strength distributions. These distributions were generated 
from fatigue data of specimens which were subjected to an alternating 
bending stress and a constant, mean shear stress. Specimens were 
subjected to different ratios of alternating to mean stress, known as 
stress ratios, at specified alternating stress levels. The specimens 
were of SAE 4340 steel and were of a grooved geometry. A grooved and 
ungrooved test specimen are shown in Figure 3.1.1 and 3.1.2. A com- 
plete description of the test program, which was sponsored by the 
National Aeronautics and Space Administration under Grant Number 03- 
002-044 at The University of Arizona, and of the procedures and 
materials is given in NASA CR-120831 (3). 

The methodology used by John Smith in developing the strength 
distribution was similar to that used by Richard E. Smith (discussed 
in Chapter II of this report) (2). Cycles to failure distributions 
were developed at specific alternating stress levels for the stress 
ratios of infinity, 3.5, 0.825, and 0.44. The lognormal distribution was 








Ungrooved Tensile Test Specimen of 4340 SAE Steel 
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found to best describe the cycles to failure data (3, p. 47). 

The decision to accept the lognormal distribution over the 
normal distribution was based on the Kolnogorov-Snirnov test 
(3, p. 44). 

Cycles to failure distributions for particular stress 
ratios at specified stress levels were plotted on S-N diagrams 
of alternating stress versus log cycles. This specified the 
mean line and +3<5" ( & - standard deviation) envelope. After 
this was established, it was possible to interpolate many cycles- 
to-failure distributions at intermediate stress levels. This 
interpolation process made it possible to calculate strength 
distributions at specific cycle life values, 

A review of the methodology used in the calculation of 
strength distributions at specific cycle life values is in ord- 
er at this time. A histogram can be located along the N 
cycle life line in such a way that the midpoints of the cells 
are at the interpolated stress levels. (See Figure 3.1.3). 

The ordinate of each of the strength histogram cells is the 
area bounded by the N cycle line and the cycle to failure dis- 
tributions which has been interpolated for that particular 
value of alternating stress. Thus if f (N/S^) is the cycle to 
failure distribution at a particular alternating stress level, 
the ordinate of the strength histogram cell at that stress 
level will be given by: 

n 

J f(N/S ± ) dn 

■ -oo 


FCN'S^ = 


(3.1.1) 
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The ordinate of each of the strength histogram cells w ill 
be : 

n 

F(N/S i ) = J f(U/S i ) dn (3.1.2) 

-<x> 

The histogram which is developed in this way is ‘the 
cumulative strength histogram of specimens failing by N 
cycles. If S is the strength variable along the H cycle life 
line, the probability density function can be developed in 
the following manner. The value of the i^ 1 cell of the 
strength probability density histogram is given by: 

f (S. ) = F(N/S. ) - F(N/S,.~1) (3.1.-3) 

11 * 

The fact that there are many interpolated cycles to failure 
distributions on the N cycle life line insures an adequate 
number of class intervals in the strength histogram. A 
normal distribution is then fitted, by statistical methods, to 
the strength probability density histogram (3, p. 55). 

3.2 Computer Method 

John Smith developed two computer programs idiich were 
used in the strength distribution calculations. The first of 
which, known as CYTOFR, calculated the cycles to failure dis- 
tribution parameters, mean and standard deviation as veil as 
the coefficients of skewness and kurtosis of the normal and 
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lognormal distributions. It also performed the Chi -Squared and 
Kolmogorov-Smirnov goodness of fit tests to determine if a normal 
or lognormal distribution fits the data best. A flow chart, 
variable definition and computer listing of this modified to 
include a sort routine is given in Appendix C. 

The program STRENG finds the normal strength distributions 
from the lognormal cycles to failure distribution parameters. Smith 
used this method because it has been found by earlier studies that 
the normal distribution adequately describes the strength data (2) 

(4). The input data of this program includes the cycles to failure 
data, two extrapolated lognormal distributions on either side of the 
experimental distributions and interpolated lognormal cycles to 
failure distributions between the actual experimental failure dis- 
tributions. The program then calculates the mean, standard deviation, 
and coefficients of skewness and kurtosis for the normal strength 
distributions. In addition a goodness of fit test, the Kolmogorov- 
Smirnov test is performed on the normal strength distributions. 

A -flow chart, variable definitiions computer listing of program 
STREX'G given in Appendix D. 
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3.3 Results 

The statistical S-N diagrams, which present the cycles to 
failure distributions for stress ratios of infinity, 3.5, 0.825 and 
C .44, are given in Figures 3.3.1 through 3.3.4. The mean, standard 
devistions, and three sigma limits of the lognormal cycles to 
failure distributions are presented in Table 3.3.1. 

John Smith presented only one S-N diagram which had the cal- 
culated strength distribution for a stress ratio of 3.5 and cycle 
life values of 10,000, 50,000, and 100,000 cycles. The normal 
parameters of the strength distribution which was placed on this 
particular S-N diagram are given in Table 3.3.2 while the S-N diagram 
appears as Figure 3.3.5. The completion of the stress ration tests 
of 0.44 allowed the calculation of strength distributions for this 
stress level. Table 3.3.3 presents the parameters of the normal 
strength distributions at three cycle life values and a stress ratio 
of 0.44 while these distributions were added to Figure 3.3.4. 

A complete table of the normal strength distribution parameters 
at all stress ratios and cycle life values was inavailable and con- 
sequently recovered by use of the program STRENG. The results 
are presented in Table 3.3.4. 
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Table 3.3.2 Strength Distribution Parameters of the Normal Distribution 
at Three Cycle Life Values For A Stress Ratio of 3.5. 

( 3 , p • 5 7 ) 


Cycles 

Parameter Estimate of 
Normal Distribution 

of 

Life N 

Mean (psi) 

Standard Dev. 
(psi) 

10,000 

106,639 

3,256 

50,000 

79,021 

3,253 

100,000 

70,972 

2,313 
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Table 3.3.3. Strength Distribution Parameters of the 
Normal Distribution at Three Cycle Life 
Values for a Stress Ratio of 0.44. 







Table 3.3.4 Parameters of Normal Strength Distribution 

At Specific Stress Ratios and Cycles of Life 



5 , 
7 , 
9 , 
10 , 
20 , 
30, 
40, 
50, 
60, 
70, 
80, 
90, 
150, 
I 200, 
1 , 000 , 


,000 
', 000 * 
',000 
1,000 
1,000 
1,000* 
',000 
',000 
>,000 
>,000 
> , 000 * 
1,000 
, 000 * 
>,000 


135,953 

128,356 

121,191 

115,926 

113,844 

101,523 

94,972 

90,693 

87,391 

84,707 

82,484 

80,628 

79,054 

72,901 

70,172 

56,184 


Standard 

Deviation 

psi 


1,629 
2,570 
' 2,628 
2,790 
2,937 
2,907 
2,801 
2,772 
2,559 
2,665 
2,671 
2,697 
2,662 
2,594 
2,539 
2,133 
1,881 
2,138 


-3 Sigma 
Limits 
psi • 


159,932 
131,533 
128,070 
119, 9S6 
112,378 
107,205 
105,441 
93,206 
87,296 
82,878 
79,378 
76,617 
74,498 
72,845 
71,437 
66,502 
64,529 
49,771 


♦Parameters used in Chapter V Finite Life Gooaman Diagrams 













Table 3.3.4 Paraneters of Normal Strength Distribution 

At Specific Stress Ratios and Cycles of Life 
(Continued). 

R = 3.5 


Cycles 


Mean Standard -3 Sigma +3 Sigma 

Strength Deviation Limits Limits 

psi psi psi • psi 


1,000 

158,162 

3,000 

132,668 

3,500* 

128,964 

5,000 

120,436 

7,000 

113,156 

9,000* 

108,523 

10,000 

77,689 

' 20,000 

94,273 

30,000 

87,100 

40,000* 

82,345 

50,000 

79,021 

60,000 

76,533 

70,000 

7 4,629 

80,000 

73,152 

90,000* 

71,967 

100,000 


150,000 

67,341 

200,000* 

64,792 

000,000 

50,533 


2,144 

3,536 

3,643 

3,797 

3,325 

3,213 

2,502 

3,611 

3,701 

3,467 

3,253 

3,049 

2,792 

2,554 

2,393 

2,327 

2,400 

2,812 


151,731 

122,061 

118,035 

109,645 

103,178 

98,883 

83,440 

75,998 

71,942 

69,261 

67,385 

66,253 

65,490 

64,787 

60,360 

57,593 

42,097 


164,592 

143,275 

139,894 

131,826 

123,133 

118,163 

105,106 

98,203 

92,747 

88,781 

85,681 

83,005 

80,814 

79,147 

74,321 

71,991 

58,969 


♦Parameters used in Chapter V Finite Life Goodman Diagrams, 
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Table 3.3.4 Parameters of Normal Strength Distribution 

at Specific Stress Ratios and Cycles of Life. 
(Continued) . 



♦Parameters used in Chapter V Finite Life Goodman Diagrams. 







1 

3 

3 

5 

7 

9 

10 

20 

30 

40 

50 

60 

70 

80 

90 

150 

200 

1,000 


,500* 

,000 

,000 

, 000 * 

,000 

,000 

,000 

1 , 000 * 

,000 

,000 

,000 

,000 

, 000 * 

,000 

, 000 * 

,000 


154,197 
130.027 
126,369 
117,923 
110,448 
105,645 
103,725 
91,498 
85,248 
80,953 
77 703 
75,179 
73,148 
71,443 
69,970 
64,263 
61,737 
40,686 


2,056 

3,651 

3,739 

3,893 

3,539 

3,308 

3,323 

3,243 

3,148 

3,208 

3,154 

3,p44 

2,974 

2,947 

2,927 

2,372 

2,133 

2,394 


148,030 

119,073 

115,151 

106,229 

99,831 

95,722 

93,756 

81,769 

75,805 

71,327 

68,239 

66,045 

64,226 

62,602 

61,189 

57,147 

55,338 

41,505 


+3 Signa 
Limits 
psi 


160,364 

140,981 

137,587 

129,618 

121,065 

115,569 

113,694 

101,227 

94,691 

96,578 

87,166 

84,312 

82,070 

80,285 

78,750 

71,379 

68,135 

55,868 


^Parameters used in Chapter V Finite Life Goodman Diaarans 
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3.4 Discussion as to Validity 

The methodology which John Smith presented was found to 
be well grounded and accurate. Disturbing to the original 
investigator was the fact that the vertical strength distri- 
butions did not completely "fill" or span the entire width of 
the statistical S-N diagram envelope. That is, the vertical 
strength distributions did not span the complete distance be- 
tween the plus three signa line and minus three sigma line of 
the cycles to failure distributions. This, however, was not 
disturbing to this investigator for the following reason. 

The histogram of the vertical strength distributions 
was developed from the cumulative failure probabilities of a 
large number of cycles to failure probability density functions. 
Along a particular cycle life line on the statistical S-N 
diagram, See Figure 3.4.1, the vertical distance from the 
intercepted plus and minus three sigma lines is of no particu- 
lar significance. 

The fact that the vertical strength distribution does 
not span the distance between point A and B can be explained 
in the following manner: Hie cycle s-to-failure distributions 

are lognormal while the strength distributions are normal. It 
must be noted that in actuality only distributions to the left 
of the cycle life line contribute percentile areas to the 
cumulative strength distribution. The two points, A and B, 
are then completely unrelated. The plus and minus three signa 
limits, which are not necessarily points A and B, of the 
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vertical strength distributions are derived iron the cumulative 
failure probabilities of the individual cycle to failure proba- 
bility density functions. The plus and minus three sigma 
limits of the vertical strength distribution are related to 
the variability of the individual probability density functions 
of the cycles to failure data at particular stress levels 
intercepted by the vertical cycle life line. 


CHAPTER IV 


GENERATION OF FINITE LIFE GOODMAN DIAGRAMS (METHOD I) 

4.1 Theory 

To date most of the work done with the modified Goodman 
Diagram and Goodman surfaces has been concerned with the 
objective of presenting the relationship of alternating and 
mean stress for infinite periods of life. For the biaxial 
stress condition, where the specimen is subjected to both 
bending and shear stress, the Goodman diagram represents 
combinations of these two stresses, in the cartesian plane, 
where the specific combinations will not cause fracture to 
occur over a period of infinite life. The need for a surface 
to describe this relation occurs where probabilistic methods 
are used in reliability calculations by the interference 
method. 

. It is the purpose of the method under discussion in 
this Chapter to break away from the traditional concepts of 
using a Goodman surface to graphically represent the relation 
between bending and shear stresses for a period of infinite 
life, and present, in the cartesian plane, the relation of 
alternating stress to mean stress where the alternating 
stress is a bending stress while the mean stress is a shear 
stress for finite periods of life. Hie advantage of such an 
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approach is that it will graphically illustrate to a designer 

that it is possible to have higher combinations of bending and 

shear stresses for a finite life design. Hence, if a part 

need only fuction for a specified finite number of hours or 

cycles, after which its failure is not detrimental to the 

success of the mission, it can be subjected to combinations 

of bending and shear considerably higher than if it had 

6 

to function in excess of 10 cycles. 

The data which was used for this method was generated by the 
Reliability Research Laboratory of the Aerospace and Mechanical 
Engineering Department under NASA Grant 03-002-044 at Die 
University of Arizona. Cycles to failure data was generated 
for rotating specimens which were subjected to an alternating 
bending and constant shear stress at specific stress ratios. 

Stress ratio is defined as the ratio of alternating bending 
stress, t.o mean normal stress from torque. 

The cycles to failure data was generated from 1967 to 
1970 for stress ratios, of infinity, (pure bending), 3.5, 

0.8, and 0.44. Die specimens which were subjected to a 
bending and shear stress had a grooved geometry. The 
material which the specimens were made of war. SAE 4340 steel. 

A comprehensive explanation of the test program, machines, 
procedures and materials used can be found in NASA CR-72839 
(3). The data generated, by this experimental effort, which was 
used to construct finite life Goodman surfaces displayed in this 
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Chapter, included the bending and shear stress in each of the over 
300 test specimens which were run to failure. The value of the 
cycles life at which the specimen failed was also recorded. 

In the original test program twelve to eighteen test specimens 
were run at each stress level for each stress ratio. The lower 
number of specimens, twelve, occurred at higher stress levels, where 
variability in the cycle life data was small. Hie larger sample 
size was used at lower stress levels where variability in the data 
suggested a larger number of specimens be used. Referring to 
Figure 4.1.1, it can.be seen that if the infinity stress ratio, lo- 
cated along the alternating stress axis, is to be bridged to a mean 
stress distribution at a stress ratio of zero by a Goodman surface, 
for the finite periods of life, there will have to be at least two 
other stress ratios which contain enough data points to form two 
strength distributions between the alternating stress distribution 
and the mean stress distribution and the mean stress distribution. 
Because of this restriction, it was necessary to screen the cycles 
to failure data to determine appropriate cycle life ranges. This 
was accoiiplished with the aid of Table 4.1.1. This table records 
the... cycles to failure data fbr specimens at particular stress ratios. 
The cycle life ranges which were determined by this method are as 
follows: 
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Table 4.1.1. Screening of Cycles-to-Failure Data to Ascertain 
Appropriate Cycle Life Ranges. 


/cM ' 


iP ~ * . x*, - * * , ** * V* ***** 

'* >x ***** ■ t, * x x x 


_ 120,000 


S 2 


*x> 

* * • 
/ 


ju. 




A K, 


x »*/" / 

*« • /X 


x V ,* x 

'** V 


1,000 


c 


49 


N " 1: 1,000 - 3,500 cycles to failure 

- 2: 20,000 - 40,000 cycles to failure 
N = 3: 60,000 - 90,000 cycles to failtnre 
N = 4: 90,000 - 200,000 cycles to failure 
N = 5: 6,000 - 9,000 cycles to failure 

It can be seen that these groups span the cycle life 
spectrum from low to relatively large cycle lifes. However, 
all five groupings are in the cycle life range to the left 
of the “knee" of the S-N diagram. It would have been desirable 
to keep the cycle life groups as small as possible, however 
the restriction of having enough data points within each 
group to form a distribution was- the governing restriction 
in this case. The alternating stress distribution, where the 
shear stress or mean stress is zero, was obtained directly 
from the cycle to failure data of pure bending, r = oo, 
specimens. After the cycle life qroups were determined the 
distribution placed on the alternating stress axis for each group 
was specified by the bending stress recorded in each of the 
.. specimens which failed within the cycle life range of that 
particular group. 

The distribution which was used on the mean stress axis 
in all cases was the ultimate strength distribution of an un- 
grooved tensile test specimen which had the same cross-sectional 
areas as the actual grooved fatigue test specimens (grooved 
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specimen results are used in Chapter V’> see Figures 3.1.1 and 
3.1.2 . The mean normal strength distribution was obtained 
directly from static tensile strength tests. In the case of 
such tensile tests, it is found that the grooved specimens 
have a higher ultimate strength. This is caused by a radial 
stress which is introduced into the specimen at the root of 
the groove. However, a grooved specimen which is subjected 
to static torque load would not experience such a radial 
stress At this time in the investigation, it would, therefore, 
not have been correct to use the grooved specimens to determine 
the strength distribution for the mean stress axis (see 
Chapters V and VII). 

After obtaining the cycle life groups the next step 
was to relate the shear stress in each of the specimens to the 
proper mean stress. This was done by using two predominant 
strength theories. These two theories are the Von-Mxses 
Hencky theory and the maximum shear stress theory. A 

complete discussion of these two theories appears in 
Chapter VIII of this report . 

The resultant stress vector, S ? , for each of the data 
points has a mean-stress of S n = ( 'XS s shear stress), 

(3, p. 3). The resultant stress vector is a combination of 
mean and alternating stresses. Tho magnitude of the resultant 


vector is; 
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Sj ■ ( S Q ) 2 + ( s m >2 

’ 1/2 

Sj= ^_(S 0 ) + ( f3t ) 2 } (4.1.2) 

In the case where the Von Hises Hencky criterion is used 
to relate the shear stress to the mean stress the resultant 
stress vector is called stress vector I. 

If the maximum shear stress theory is used to relate the 
mean stress to the torsional stress, the governing relation 
will be (13, p. 2): 

S = 2^ 
m 

The maximum shear stress theory predicts that yielding 
will occur when the maximum shear stress is equal to the shear 
stress corresponding to the shear stress produced in a simple 
tension test for yield strength (15, p. 152). The Mohr circle 
predicts that yielding will begin when (15, p, 152): 

^max = V 2 (4.1.3) 

Hence, the mean stress is given by S B = 2 

Although both these failure theories are based on yielding 
as the failure criterion much experimental data indicates that 
they apply as well when fracture is the failure criterion as 
represented by the static ultimate strength. 


stress 

(4.1.4) 

(4.1.5) 

If the maximum shear stress theory is used to relate the shear 
stress to the mean stress the resultant stress vector is referred 
to as stress vector II. 

The resultant stress vector must be described by both a 
magnitude and a direction. A stress ratio was previously 
defined as the ratio of alternating stress to mean stress. 

It can be seen that the stress ratio will vary with the 
strength theory which relates the shear stress to mean stress. 
Hence, for the stress vector where the Von Mises Henchy theory 
was used the stress ratio becomes: 

r 1 = Sa / 'T3'£ (4.1.6) 


For the maximum shear stress theory the resultant 
vector magnitude will be giv;n by; 


=7 « < S a > 2 + < S B > 2 


J = [is/ ♦ (2'C) 2 J 


1/2 


For the second stress vector in which the maximum shear stress 
theory is used to relate shear stress to mean stress: 
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For both stress ratio cases the data had variability, 
described by the standard deviation of r (CT_ .) , which was 
quite small. The mean values of r and the corresponding 
standard deviation of each of the cycle life groups is given 
in Table 4.1.2. Because the variability of r was small in 
each case the resultant stress vector was assumed to lie on 
♦he stress ratio axis, defined by the average value of r. 
Referring to Figure 4.1.1 it can be seen that the angle G> , 
along which the resultant stress vector is oriented is 
given by: 

—1 

© = tan (S /S ) (4.1.8) 

Q W 

The orientation of the stress vector is thus specified 

by the mean stress ratio, HT , for each cycle life group. 

For each data point in a cycle life group, it was possible to 

use the PDP-S Computer to perform the calculation to obtain 

values for S r . After the S r values specified by Equations 

(4.1.2) and (4.1.5) were obtained for each of the data points, 

it was possible to compute a mean S r and a standard deviation 

(ct' c ) which then specified a strength distribution. Ihese 
* 

were then plotted along the mean r-axis. After the values of ST 

r 

and 0"g were obtained the Kolraogorov-Snirnov test was used 
to determine if the normal or the lognormal frequency functions 
could be accepted as representing the strength distributions. 


A 
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Table 4.1.2 Mean and Standard Deviation of Stress Ratios. 


Cycle Life Range 

Mean 

T 

Standard Deviation 

cr~ 

r 

. N b 1,000-3,500 



Stress Vector I 

3.514 

0.1200 

Stress Vector II 

3.043 

0.1306 

N = 6,000-9,000 



Stress Vector I 

3.498 

0.1490 

1 

0.875 

0.0290 

Stress Vector II 

3.029 

0.1296 


0.757 

0.0254 


N = 20,000-40,000 
(. Stress Vector I 

Stress Vector II 

H = 60,000-90,000 
Stress Vector I 

Stress Vector II 


3.508 

0.1430 

0.729 

0.0890 

3.043 

0.1239 

0.631 

0.0188 


3.422 

0.1720 

0.872 

0.0910 

2.963 

0.1492 

0.755 

0.0171 


H = 90,000-200,000 

Stress Vector I 0.804 

0.439 

Stress Vector II 0.696 

0.382 


0.0556 

0.0710 

0.0482 

0.0110 



\ - 
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4-.' 2 Computer Method 

Three computer programs of varying complcx5.ty were used 
to determine the strength distributions discussed in this 
chapter. The first of these, a PDP-8 computer program referred 
to as BAR I, Tised an input of bending stress and shear stress 
to obtain a value of S r when the Von 1-fises Hencky strength 
theory was used. For a shear stress (and alternating stress) 
it performs the following calculations. It first calculates 
a mean stress from a shear stress: 

S = fs't (4.2.1) 

Id 

and then performs the operations required by Equation 
4.1.2: 

1/2 

s r = (4.2.2) 

The second computer program referred to as BAR II which 
is also a PDP-8 program, performs the calculations required 
by the maximum shear theory to relate the shear stress to mean 
stress: 

*, » (4.2. 3) 

BAR II then performs the calculations required by Equation 
4.1.5: 


s? “ f(s 0 ) 2 + (a*) 2 } 


1/2 


(4.2.4) 


A flow chart, variable definitions and computer listing of the 
PDP-8 programs BAR I and BAR II is given in Appendix E. 

After the individual values for the resultant stress 
vector Magnitudes for each specimen were computed the data 
was then submitted to the CDC-6400 program ASMERG which 
made the foilwing calculations. The mean and standard devia- 
tion of for each cycle life group ‘and stress ratio was 
calculated. These values are denoted by S r and o . In 
addition the Kolmogorov -Smirnov goodness of fit test, for 
a normal and lognormal distribution, was conducted on each of 
the strength distributions. The coefficients of skewness and 
kurtosis were determined for eaqh strength .distribution. A 
flow chart, variable definitions and computer listing of program CYTOFR 

is given in Appendix C. 

4.3 Results 

For the two strength theories considered, Von Mises 
Hencky and maximum shear stress, the results of the calcula- 
tions to determine the resultant stress vector magnitudes are 
given in Table 4.3.1. In addition these tables present the 
stress ratio mean for each cycle life group as well as the 
calculated value for the mean stress for each strength theory 
considered. Tables 4.3.2 and 4.3.3 present the mean and 
standard deviation of the resultant stress vector magnitude 
for each stress ratio by cycle life group. Table 4.3.4 


presents the ultimate strength distribution of an ungrooved 
test specimen. Table 4.3.5 presents the Kolnogorov-Srairnov 
/D/ values for each cycle life group for both the Von 

EtQX 

Hises-Henchy strength theory as well as the maximum shear 
stress theory. 

Based on the Kolmogorov- Smirnov test, in.no case can the 
proposed distributions, normal or lognormal, for the result- 
ant stress vector be rejected at the ninety percent confidence 
level. In ten out of fifteen resultant stress vector distri- 
butions the lognormal distribution had smaller /D/_„, values • 

For the normal and lognormal the / D/ difference was quite small 

max 

the difference occuring usually in the third decimal place. 
Because of this slight difference, it cannot be said that 
either the normal or lognormal distributions .fit the data bet- 
ter, rather that both the normal and lognormal distributions 
fit the strength data equally well. 

The third and fourth statistical moments, coefficients 
of skewness and kurtosis, give little insight into the nature 
of the underlying distribution. The "knowledge of the third 
moment gives almost no clue as to the shape of the distribu- 
tion" (5, p. 109). 



Table 4.3.1 C/cles to Failure .Data, Bending and Shear Stress Data and the Resultant 
Stress Vector Magnitudes S r (I) and S r (II) for the Von Mises Hencky 
and the Haxiraum Shear Stress Theory. 
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Table 4.3.1 Cycles to Failure Data, Bending and Shear Stress Data and the Resultant 
Stress Vector Hangitudes S r (I) and S r (II) for the von Mises Hencky 
and the Maximum Shear Stress Theory. 
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Table 4.3.1 Cycles to Failure Data, Bending and Shear Stress Data and the Resultant 
Stress Vector Magnitudes S r (I) and (II) for the von Mises Hencky 
arjd the Maxirvoei Shear Stress Theory. 
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Table 4.3.1 Cycles -to Failure Data, Bending and Shear Stress Data and the Resultan- 
Stress Vector Magnitudes (I) and S (II) for the von Kises Hcncky 
and the Maximum Shear Stress Theory. 
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0.758 256 129,228 62,835 47,856 82,839.0 104;014 95,712 114,495 
0.785 257 101,193 62,479 45,945 79,579.1 101,175 91,390 111,13.9 
0.723 258 167,735 61,711 43,956 84,794.3 104,373 97,912 115,237 


Table 4.3.1 Cycles to Failure Data, Bending and Shear Stress Data and the Resultant 
Stress Vector Magnitudes S r (I) and S r (II) for the von Mines Hencfcy 
and the Maximum Shear Stress Theory. 
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Table 4.3.2 Mean/ Standard Deviations and + 3 Limits of 
Resultant S .ress Vector for the" von Mises- 
Henclcy Theory. 


Life Range - Cycles 

Mean 

psi 

Standard 

Deviation 

psi 

+3®’ 

psi 

-3 cr 
psi 

11=1,000-3/ 500 

r = cd 

144,327 

1,538 

148,934 

139,710 

r e 3.514 

156/651 

3,683 

167,700 

145,602 

r s 0.441 

150,404 

3,261 

170,007 

150,801 

K=6, 000-9,000 

r = oo 

113,862 

'954 

116,697 

111,027 

r = 3.498 

119,228 

1,999 • 

125,225 

113,231 

r = 0.875 

168,382 

3,170 • 

177,892 

158,872 

N=20, 000-40, 000 

r = cd 

98,247 

2,683 

106,300 

90,196 

r = 3.507 

86/301 

1,062 

89,487 

83,115 

r = 0.729 

146,952 

3,883 

158,701 

135,203 

H=60, 000-90, 000 

r = co 

81,282 

993 

84,821 

78,303 

x = 3.422 

78,211 

2,550 

35,861 

70,561 

r = 0.872 

117,097 

3,900 

123,797 

105,397 

R=90, 000-200, 000 
r = oo 

74,336 

3,556 

85,004 

.'63,668 

r = 0.804 

104,074 

8,517 

112,591 

95,557 

x = 0.439 

148,137 

2,325 

155,112 

141,162 
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The Goodman surfaces and dio.rjrc.ns generated by this method 
are presented in Figure 4.3,1 through Figure 4.3.10 for the various 
cycle life groups and strength theories. 

4 . 4 Discussion as to Validity 

The Goodman surfaces generated illustrate that for 
cycle life ranges which are well below infinite life the 
allowable combinations of bending end shear stress are much 
larger in magnitude than combinations of the same biaxial 
stress which could be sustained by a specimen for an infinite 
life range (Figure 4.4,1 is a Goodman Surface for a period 
of infinite life). As the cycle life groups mean value in- 
creases the corresponding combination -of bending and- shear 
stress continues to decrease toward the value which is presented 
in the infinite life diagram. A comparison of this surface to 
the surfaces generated for finite life periods, see Figure 4.3.9 

illustrates the above conclusion. 

The variability in the standard deviations of the resultant 
stress vector for stress ratios other than infinity and zero 
should be examined closely. This variability 

is larger in many cases than the variability of the bending 
stress and ultimate strength distributions used for the alter- 
nating and mean stress axis. In reviewing the methodology 
which was discussed^ one finds that data used for the alternating 
stress distribution and the ultimate strength distribution 
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were taken from one series of tests. However, the data which was 
used to generate the strength distributions for ratios of 3.5, 
0.825, and 0.44 were taken from a series of tests, each with their 
own variability. From the algebra of normal functions we sec that 
when two standards deviations are added the resultant standard 
deviation becomes (16, p. Ill): 

I 2 2 ~ ~ 

a =v/ 0 + o +2po 

x+y V x y x y 

where p is the correlating coefficient. If p = 0, assuming 
independence (16, p. Ill): 

x r~2 2 

a = V o +0 

x+y V x y 


(4. '.2) 


N=1 1000-3500 CTF 

Stress Vector T 


S 



Fig. 4,3.1 finite Life Goodrian Diagram and Surface 



Fig* 4.3.2 Finite Life Goodman Diagram and Surface. 



Fig. 4.3,'3 'FinlteLife Goodman' Diagram arid Surface. 
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3.6 Finite Life Goodman Diagram and Surface. 
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Fig. 4.3.7 Finite Life Goodman Diagram and Surface 



.3.8 Finite Life Goodman Diagram and Surface. 
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Table 4.3.3 Mean, Standard Deviation and + 3 Limits of 
Resultant Stress Vector for the Maximum Shear 
Stress Equation. 


Life Range - Cycles 

Mean 

psi 

Standard 

Deviation 

psi 

+3 o' 
psi 

-3c- 

N=1 .000-3 .500 

r . ~ 3.043 

158,600 

3,641 

169,523 

147,677 

r = 0.381 

181,411 

3,704 

192,523 

170,299 

N=6, 000-9, GOO 

r = 3.029 

120.787 

2,007 

126,808 

114,766 

r = 0.757 

183,578 

3,812 

195,014 

172,142 

H=20,000-40>000 
r = 3.043 

87,375 

1,077 

90,606 

84,144 

r = 0.631 

162,152 

4,178 

174,686 

149,518 

N~60,000-90,000 

r = 2.960 

79,236 

2,586 

89,994 

71,478 

r = 0.755 • 

127,704 

4,405 

140,916 

114,492 

N=90, 000-200, 000 
r = 0.696 

114,126 

2,896 

122,814 

105,438 

r = 0.332 

167,554 

2,742 

175,780 

159,328 


Table 4.3.4 Ultimate Strength Distribution, Mean, Standard 

Deviation, and + 3. Limits for Moan Stress Axis 
for Ungrooved Specimens . 


Mean 

Standard +3o 

Deviation 

-3 a 

psi 

psi psi 

psi 


r = 0 


176,000 


2,500 


185,500 


170,500 
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Table 4.3.4 Ultimate Strength Distribution, Mean, Standard 
Deviation, and +3 Limits for Kean Stress 
Axis for Ungrooved Specimens. 


Mean Standard +3<r -3 a* 

Deviation 

psi psi psi psi 

2,500 


1=0 


173,000 


185,500 170,500 


Table 4.3.5 Kolnocorov-Snirnov /D/ Values for Stress Vectors I and II 
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CHAPTER V 


■ 'GENERATION op FINITE LIFE GOODMAN DIAGRAM (METHOD II) 

5.1 Theory 

It is possible to place strength distributions which 
have been developed iron cycles to failure distributions on 
statistical Goodman diagrams, as well as S-N diagrams. In 
Chapter III, the completed results of strength distributions 
developed by John Smith for the statistical S-N diagram 
were presented. 

It is possible to place the strength distributions 
developed by the technique discussed in Chapter III on finite 
life. Goodman diagrams . It is first noted that the vertical 
strength distributions developed in Chapter III are for a 
specific cycle life. This finite life Goodman diagram is then 
for this cycle life. 

Figure 5.1.1 compares a strength distribution placed on 
a finite life Goodman diagram to that of the same strength dis- 
tribution placed on a statistical S-H diagram. Tne ordinate 
axis of both the Goodman diagram and the statistical S-H dia- 
gram is the alternating stress level. Because of this the 
vcilue of the mean of the alternating strength distribution 
when transformed from an S-N diagram to a finite Goodman dia- 
gram will not change. The stress ratio and alternating stress 


alternating stress 



Ni =1000 cycles. 


N=1000 cycles 


TRANSFORMED STRENGTH 
DISTRIBUTION 


MEAN STRESS 

Fig 5.1.1 Comparison of Vertical Strength Distribution to A Strength Distribution 
Transformed to the Finite Life Goodman Diagram. 


£ 
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level, which are specified by the statistical S-N diagram, 
define, the location of the mean of the strength distribution 
which is transformed to the finite life Goodman diagram. 

The three sigma limits of the transformed strength dis- 
tribution arc not the same as that of the vortical strength 
distribution. Figure 5.1.2 depicts a strength distribution 
placed on the familiar coordinate axis of the Goodman diagram. 
The three sigma limits of the strength distribution placed on 
the stress ratio axis can easily be derived from the three 
sigma limits of the vertical strength distribution in the 
following manner. 

Equate the upper and lower three sigma limits of the 
vertical strength distribution to S. The stress ratio, r, is 
equal to the alternating stress d.ivided by the mean stress 
which referring to Figure 5.1.2 is equal to tan Q . 


r - S /S = tan© (5.1.1) 

an v ' 


-1 

© = tan (r) (5.1.2) 


The upper and lower three sigma limits of the transformed 
strength distribution are equal to S'. From Figure 5.1.2 it 
can be seen that: 


sin © = S /S' 

a 


(5.1.3) 



th Distribution s 
s to the Stress Ratio Axis. 


i 
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Solving equation 5.1.3 for S' and substituting equation 5.1.2 
yields : 

S' = S^/sin 0 (5.1.4) 

S' = S /sin { tan^Cr)) (5.1.5) 

3 - 

If. the value of r is known and the original value of the 
upper and lower three sigma limits are specified by the 
vertical strength distribution then the upper and lower 
three sigma limits of the transformed strength distribution 
can be placed along the stress ratio axis and are specified 
by equation 5.1.5, 

Five finite life Goodman diagrams were developed by the 
transformation of vertical strength distribution to the 
stress ratio axis of the Goodman diagram. The finite life 
Goodman diagrams were developed for cycle lifes of 3 / 500 / 9,000, 
40,000, 90,000 and 200,000 .cycles. Strength distributions 
placed on these finite life Goodman diagrams ware at stress 
ratios of infinity, 3.50, 0.825, and 0.44. As discussed in 
Chapter VXX of this report the strength distribution placed 
on the mean stress axis was taken to be that of the ultimate 
strength distribution of the grooved test specimen. 

A PDP-8 computer program, ROTO, was developed to perform the 
calculations required by the transformation of the upper and lower 
three sigma limits as discissed in section 5.1. A flow chart, 


variable definitions and computer listing of the ROTO program is given 
in Appendix E. 

5.2 Results 

The finite life Goodman diagrams developed by this method showing 
the strength surfaces for discrete cycles to failure are given in 
Figures 5.3.1 through 5.3.5. In Table 5.3.1 the alternating stress level 
of the mean cf the strength distributions for each stress ratio as well 
as the transformed upper and lower three sigma limits of the strength 
distributions are given. The original parameters of the vertical strength 
distributions appear in Table 3.1.3 as starred quantities, as they were 
developed by this investigator to correspond with the cycle life values 
of Chapter IV. The ultimate strength distribution of the grooved and 
ungrooved test specimen are compared’ in Table 5.3.2. 

5.3 Discussion as to Validity 

The principle of transforming vertical strength distributions to 
the Goodman diagram as explained in Section 5.1 is a Straight forward 
procedure. The transformation is simply projection of a known distri- 
bution to a different plane which in this particular case is the stress 
ratio axis of the finite life Goodman diagram. 

The vertical strength distributions have associated with them a 
cycle life value. It is possible because of this fact to develop finite 
life Goodman diagrams and surfaces from the valid vertical strength dis- 
tributions. This method is consistent, in that the vertical strength 
distributions are all developed in the same manner, and hence, there is 
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no problem of differing variabilities caused by inconsistencies in the 
procedure of developing strength distributions at various stress levels. 
The elimination of the inconsistencies in procedure in developing 
strength distributions at various stress ratios is the principle ad- 
vantage of rotating vertical strength distributions to the finite life 


Goodman diagram. 







>I«K 
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Table 5.3.1 Alternating Stress Level of the Kean and Standard 
Deviation, and + Three Sigraa Limits of Strength 
Distributions Placed on Stress Ratio Axes of the 
Finite Life Goodman Diagram. 




S 

a 

°s 

a 

-3cr 

+3 <r 

N 

- 3,500 
oo 

135,953 

2,628 

. 128,070 

143,837 


3.5 

128,964 

3,643 

118,035 

139,894 


0.825 

126,369 

3,739 

115,151 

137,587 


0.44 

105,107 

1,260 

101,327 

108,887 

N 

= 9,000 
oo 

115,926 

2,907 

107,205 

124,647 


3.5 

108,523 

3,213 

98,883 

118,163 


0.825 

105,645 

3,308 

95,722 

115,569 


0.44 

92,390 

2,002 

• 86,384 

98,396 

N 

= 40,000 



* 



oo 

90,693 

2,605 

82,878 

98,507 


3.5 

82,345 

3,457 

71,942 

92,747 


0.825 

80,953 

3,208 

*1,327 

90,578 


0.44 

72,260 

3,269 

62,452 

82,068 

N 

= 90,000 
oo 

79,054 

2,539 

71,437 

86,671 


3.5 

71,967 

2,393 

64,787 

79,147 


0.825 

69,970 

2,927 

61,189 

78,750 

N 

= 200,000 
oo 

70,172 

1,881 

64,529 

75,815 


3.5 

64,792 

2,400 

57,593 

71,991 


0.825 

61,737 

2,133 

55,338 

68,135 


0.44 

57,521 





Table 5.3.2 Comparison of Grooved 
Strength Distribution, 

and Ungrooved Specimen’s Ultimate 

Mean 

. Standard Deviation 


Grooved 255,300 psi 2,720 psi 


178,500 psi 2,500 psi 


Ungrooved 


Table 5,3.2 


Grooved 

Ungrooved 


Comparison of Grooved and Ungrooved Specimen's 
Ultimate Strength Distribution. 


Standard DeviatL on 
2,720 psi 


178,500 psi 2,500 psi 


Mean 

255,300 psi 



CHAPTER VI 

EVALUATION OP GENERATION OP FINITE 
LIFE GOODMAN DIAGRAMS ' 



6.1 Evaluation of Previous Technicpies in Developing Goodman 

Diagrams and Surfaces 

The methodology which was presented in Chapters IV 
and V for the transformation of cycles-to-failure data to 
strength distributions and the resulting Goodman surfaces and 
diagrams will now evaluated. It appears that the method of 
rotating vertical strength distrinutions and placing them on the 
stress ratio axis of a Goodman diagram is the most uniform method 
of creating the Goodman surface. As discussed in Section 5.4 all 
distributions are created in the same manner which eliminates the 
problems of variability discussed in Section 4.4. 

In viewing the general chape of both the finite life Good- 
man diagrams developed in Chapter IV and V the following con- 
clusion can be dra^m. By design each of the cycle life groups 
in both Chapters IV cmd V are similar . This was done for 
purposes of comparison. Die curves shift progressively lower 
as cycle life value increase. The curves developed in both 



chapters are consistent in this respect. As cycle life design 
values decrease larger combinations of bonding and shear stress 
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The two sets of finite life Goodman diagrams do have 
some inconsistencies, . The ultimate strength distribution of 
the ungrooved test specimen was used as the static strength 
distribution in Chapter IV while the ultimate strength dis- 
tribution of the grooved specimen was used in Chapter V» Both 
distributions were placed along the mean stress axis of the 
Goodman diagram. An explanation of these facts is in order 
at this time. 

The method used in Chapter IV to generate the finite life 
Goodman diagrams and surfaces was completed in early February 
of 1971. At that time, there was doubt as to the geometry of 
the specimen to be used, consequently the conservative ulti- 
mate strength distribution of the ungrooved specimen was 
chosen. Although this doubt did exist the diagrams developed 
in Chapter IV did not seem to support the use of the higher 
ultimate strength of the grooved geometry specimens. 

The finite life Goodman diagrams which were developed 
in Chapter V, in late March and early April of 1971, illustrated 
the critical importance of the geometry of specimen from 
which the static ultimate strength distribution was derived. 

All but the finite life Goodman diagram developed at 200,000 
cycles showed the mean stress values a.t a stress ratio of 0.44 
close to or obove the ultimate strength distribution of the 
ungrooved test specimen. This fact, although quite alarming at 
f irst, made an investigation of the correct geometry of paramount 
importance and considerable effort, as exhibited in Chapter VII, 


vas made "to resolve the question concerning the proper ultimate 
strength distribution to be placed along the mean stress axis* 

6 , 2 Recommendations 

The weaknesses and advantages of the two methods of 
generating Goodman surfaces and diagrams have been discussed. 

It' is the opinion of this investigator that the method devel- 
oped in Chapter V, that of forming vertical strength distribu- 
tions from cycles to failure data and then transforming these 
to the stress ratio axis of o. Goodman diagram, be considered 
as the appropriate method to be used when finite life Goodman 
diagrams are to be created from a large scale cycle to failure 
fatigue test program. Subsequently the results of Chapter V 
were weighted to a greater degree than those of Chapter IV 
in discussions dealing with the choice of geometry of test 
specimen used for the ultimate strength distribution in Chapter 
VII. This was also true in Chapters X and XI where the best 
empirical and theoretical math model of the Goodman diagram 
were sought. 

In conclusion, the transformation of vertical strength 
distributions to the Goodman diagram is considered as the ap- 
propriate methodology to be used when large amounts of cycles- 
to-failure data are available. Chapter XIII considers a more 
efficient plan whore it may be possible to generate Goodman 
surfaces with a minimum amount of actual fatigue testing. 


CHAPTER VII 


THE STATIC STRENGTH DISTRIBUTION TO BE PLACED ON THE MEAN 
STRESS AXIS OF FINITE LIFE GOODMAN DIAGRAMS 

7.1 Introduction to Static Strength Distribution 

In the previous chapters of this report, we discussed and 
developed several experimental methods of presenting fatigue 
data generated by relatively expensive and time consuming 
fatigue test programs. If however, this same Information could 
be extracted from uniaxial static tests of materials the sav- 
ings in time and effort would be tremendous. Chapters VIII and IX 
present several empirical and theoretical models of the Goodman 
diagram. 

The development of relations, which take advantage of 
information obtained from static tests, to model fatigue data 
require that the following topics be investigated. For both 
the grooved and ungrooved ge \ ;ry specimens which have been 
investigated, what are the equations which are used to 
determine the tensile yield, ultimate -and breaking strengths? 
Secondly, is it possible to specify a theoretical strength 
distribution, such as the Gaussian normal or lognormal distri- 
bution, to each of these quantities? Of major concern in the 
development of the Goodman diagram for the grooved geometry 
specimen, subjected to the combined stress condition of alter- 
nating bending and constant shear stresses, is the determination 
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of the strength attribution to be used along the neon stress 
axis. It must be determined what strength parameter should be 
used along the axis and if the diagram is to model the behavior 
of grooved specimens, on which geometry specimen, grooved or 
ungrooved, should this strength parameter be based. 

7 ; 2 Calculation of Yield / Ultimate and Breaking Strengths 

The calculation of static tensile strengths, is well 
documented throughout the literature. The yield strength 5.s 
calculated by dividing the force initiating the yield by the 
cross-sectional- area of the specimen. This cross-sectional 
area is based on the original diameter of the test specimen 
(6, p. 4). In the soft, ductile steels the yield strength is 
clearly marked by a yield point as shown in Figure 7.2.1a. 

In other materials where the yield point is less obvious, see 
Figure 7.2.1b, common paractice defines the yield load as the force 
which is required to give a 0.2 percent plastic offset (6,* p. 45). 

Ihe ultimate tensile strength is defined as the maximum 
load sustained by a tensile test specimen divided by the 
"original" cross-sectional crea. However, this calculation 
yields a parameter which is inaccurate and artificial, Ihe 
load and area on which this parameter is based do not occur 
simultaneously. For most ductile materials the maximum load 
occurs after appreciable elongation which is obviously accomp- 
anied by a reduction in area. The ultimate tensile strength, 
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Figure 7.2,1 a) Ductile Material With Clearly Defined Yield 

Point. b) Ductile Material Without Marked Yield 
Point and Comparison of True Stress-Strain Curve 
versus Engineering Stress-strain Curve. (6, p. 5) 



as calculated above, is, however, the most commonly sited 
parameter of material strength (7, p. 152). 

By definition, both the ultimate and breaking strengths 
are based on the original cross-sectional area of the test 
specimen. Hence, the breaking strength is calculated by divid- 
ing the load at fracture by the original cross-sectional area. 
Figure 7.2,1a indicates that the breaking strength of a duo- 
tile material may be less than the ultiraate tensile strength 
of the material. As necking and elongation occur in the test 
specimen the stress in the specimen continuously decreases. 

A more realistic measure of material strength is the tensile 
fracture stress. The fracture stress is determined by dividing 
the load just prior to fracture by the area measured just after 
fracture. Although the load decreases after the ultimate ten- 
sile stress is reached, the cross-sectional area decreases more 
rapidly which results in an increasing "true stress." Because 
of this the fracture stress is equal to or greater than the 
ultimate tensile stress (7, p, 154). 

Based on the above discussion, it can be concluded that 
for the grooved and ungrooved specimens the ultimate and yield 
strength calculations should be based on the original cross- 
sectional area. Hie calculation of the breaking strength should 
be based on the reduced dianeter measured after fracture. 
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7 . 3 Theoretical Strength Distributions of the Strength Parameters 

The calculations of the strength parameters yield, 
ultimate and breaking strengths have been specified, The next 
point to be investigated is to determine if these strength 
parameters exhibit a known theoretical distribution. An 
effort to determine if the normal or lognormal distribution 
v/as favored for the three strength parameters previously 

discussed was initiated. This included gathering tensile test 

„ 

data from the research effort carried out. tinder NASA Grant 
No. 03-002-044 at The University of Arizona.' Tensile test data 
from Fhase I and Phase IX of the program was gathered for 
tensile yield, ultimate and fracture strength of both the 
grooved and ungrooved test specimens. This data vros statis- 
tically reduced by the computer program CYTOFR* v/hich has the 
ability of performing the following statistical operations; 
mean, standard deviation, coefficients of skevmess and kurtosis, 
the /D/ ttaX value for the Kolmogorov-Smirnov goodness of fit 
test and the total Chi-Sauared value, V , , for the Chi- 
Squared goodness of fit test. The AS^ERG program performs these 
statistical calculations -For both the normal and the lognormal dis- 
tributions. 

Efforts to specify either the normal or lognormal distribution 
as favoring the yield, ultimate and breaking strengths of both the 
grooved and ungrooved specimens were not totally successful. The 
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Kolmogorov-Smirnov and Chi-Squared goodness of fit tests have the 
ability of only rejecting the normal and lognormal distributions as 
properly representing the data in question. Table 7.3.1 presents 
the results of the Kolmogorov-Smirnov and Chi -Squared goodness of 
fit test results. The mean and standard deviations of the strength 
parameters are given in Table 7.3.2. The breaking strength data 
in all but the Phase I specimens was found to be rejected as being 
either normally or lognormally distributed. 

Previously the Kolmogorov-Smirnov test has been used as a basis 
in determining whether a normal or a lognormal distribution fits 
the data best based on which of these distributions had the smaller 
/D/ ax value ( 3 > P- 47). However, because of the small difference 
in the maximum value of /D / for both the normal and the lognormal 
distributions, the difference occuring in the second decimal place, 
it can not be concluded whether the normal or the lognormal dis- 
tribution gives a better fit to the experimental data. It can, 
however, be said that neither the normal or the lognormal dis- 
tribution can be rejected as distributions representing the yield, 
and ultimate strength of both the grooved and ungrooved specimens. 
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Table 7.3.1 Kolmogorov-Snirnov and Chi-Squared lest Results 
for Grooved and Ungrooved Test Specimens 


Parameter 

K ^crit V rcrit 
n-r-1 

Normal 

D 

LogNor 

D 

Normal 

X 2 

LogNor 

X 2 

UNGROOVED 







Yield 







Phase I 

10 .3686, 

- 

• 2135 

.2104 


•• 

1970 

35 .2018, 

42.3 

.1093 

.1040 



1971 

35 .2018, 

42.3 

.1173 

.1187 



Ultimate 







Phase I 

10 .3586, 

•• 

.2006 

.1979 

•» 

•* 

1970 

34,. 2047, 

41.4 

.1151 

.1097 

8.099 

8.316 

1971 

35 .2018, 

42.3 

21239 

.1255 

10.172 

10.146 

Breaking 







Phase I 

10 .3686, 


.1493 

,1461 

•• 


1970 

33 .2077, 

40,3 

.2632 

.2710 

5.919 

9.139 

1971 

33 .2077, 

40.3 

.2377 

.2277 

62.579 

48.233 

GROOVED 







Ultimate 







Phase I 

10 .3686, 


.1439 

,1420 


«* 

1970 

33 .2077, 

40.3 

.0067 

.0851 

.788 

.809 

1971 

32 .2108, 

39.1 

.0907 

.0921 . 

2.222 

2.833 

Breaking 



• 




Phase I 

10 .3686, 

- 

.1229 

.1238 


•* 

1970 * 




M 

m 

m 

1971 * 



' 


m 

m 


X test could not be ‘applied due to insuff^cent data points. 


* 


Breaking load not measured with sufficient accuracy to calculate 
strength distribution. 
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Table 7.3.2 

Mean Values and Standard Deviations of Tensile 
Strength Distributions (psi) 


Mean 

Standard Deviation 

UKGEOOVED 



Yield 



Fhase I 

171,150 

2,779 

1970 

158,285 

5,840 

• 1971 

155,505 

1,765 

Ultimate 

Phase I 

177,850 

2,582 

1970 

167,044 

5,273 

1971 

165,103 

1,521 

Breaking 

Phase I 

254,800 

4,391 

1970 

255,904 

12,964 

1971 

260,921 

8,247 

GROOVED 



Ultimate 

Phase I 

255,300 

2,720 

1970 

254,380 

2,260 

1971 

269,137 

2,832 

Bre aiding 
Phase I 

303,950 

3,122 

1970 

m 


1971 

m 
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A survey of the literature clearly indicated that strength data 
is usually assumed tc be distributed normally. The NERVA Project 
Report states that, "Usually the data (strength data) will be assumed 
to be normally distributed, however, lognormal and Weibull distri- 
butions are acceptable and can be used much the same as normal data." 

(8, p. 9). Hence, it can be concluded that the normal distribution is 
acceptable as the theoretical distribution considering the limited ex- 
perimental evidence available and the fact that there is an adequate 
amount of documentation in the literature to support this conclusion. 

7.4. Mean Stress Axis Strength Parameter 

The determination of the strength parameter, and consequently, the 
strength distribution to be used in any fatigue data model is dependent 
upon the model used and the definition of the failure mode. If it has 
been determined that yielding is detrimental to the proper functioning 
of the specimen then the distribution of yield strength should be used. 
If only the fracture is of concern then the ultimate strength dis- 
tribution should be used (10, p. 4). Chapter X describes the Goodman 
line in detail. The Goodman line connects the endurance strength to the 
ultimate tensile strength. The failure criterion of the cycles to 
failure data presented in this report has been fracture. The ultimate 
strength distribution is concluded to be the proper distribution to be 
placed along the mean stress axis if the failure criterion is fracture 
and the fatigue model is the Goodman diagram. 
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7.5. Least Squares Estimate of the Ultimate Tensile Strength 

The final topic which remains to be discussed concerns the 
geometry of the tensile test specimen whose ultimate strength dis- 
tribution will be placed on the mean stress axis of the Goodman diagram. 
There exists the possibility of using the grooved or the ungrooved 
ultimate strength distribution of a tensile test specimen of equal 
cross-sectional area. 

Initially one might conclude that if the Goodman diagram is to 
model the behavior of a grooved specimen then the grooved ultimate 
strength distribution should be chosen. Mr. Carl S. Osgood, author of 
Fatigue Design , in response to a letter which solicated his opinion 
on this subject stated, "I believe it would be rather meaningless to 
try for a distribution on the axis for both types of specimens." 
Robert C. Juvinall, author of Stress Strain and Strength , in a reply 
to the same question suggests that the static ultimate strength dis- 
tribution" ... should pertain to the same notched (grooved) specimens 
as the S a ~S curve itself." Juvinal! concludes that the proper static 
strength distribution to use is that of the ultimate strength distri- 
bution aling the S^ axis as discussed in Section 7.4 of this report. 

There are, however, logical and well presented arguments 
supporting the use of the ungrooved specimen's ultimate strength dis- 
tribution. The mean stress axis is really at a stress ratio of zero, 
as the stress ratio of is zero at this point. Consequently, the 

alternating stress, S^, must be zero. In the case of combined bending 
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and shear this means that the mean stress axis is the axis of pure shear. 
From tensile tests, it has been found that the grooved geometry test 
specimen has a much higher ultimate tensile strength. Reviewing the dis- 
cussion found in Chapter IV of this report, this higher strength is 
caused by a radial stress which is introduced into the specimen at the 
root of the groove. However, a grooved specimen which is subjected to 
a static torque load would not experience such a radial stress and would 
fail at a torsional load equal to the load which causes an ungrooved 
specimen to fail. The mean stress axis in the particular case investi- 
gated can be thought of as representing the failure mode where pure 
shear is the cause of failure. Based on the above argument, it has 
been suggested that the ungrooved ultimate strength distribution be 
used on the mean stress axis of Goodman diagrams (3, p. 71). 

In attempting to resolve these two differing opinions this in- 
vestigator turned to an analytic evaluation of the problem. There are 
available two sets of Goodman diagrams, presented in Chapter IV and V 
of this report, on which to base such an analytic solution. Such an 
analytic solution was desired, as smooth curves can be drawn connecting 
the Goodman diagram data to both the grooved and ungrooved distribution 
of ultimate strength. 

The technique used was based on the method of least squares. The 
conventional method of least squares, however, was not considered as 
being appropriate or even workable for the problem under consideration. 

A conventional least squares analysis requires that the equation of the 
expected line be completely specified. Using this knowledge the method 
of least squares will fit the ’’best" polynomial to the data. 
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Before proceeding further, let us pause and review the problem 
and the available information. We have developed Goodman diagrams for 
varying cycle lives including that of infinite life, which have mean 
values of the strength distributions specified at stress ratios of 
infinity, 3.5, 0.825, and 0,44. In addition, there are also several 
theoretical equations which are know to model the curve which should 
be drawn between the mean of these strength distributions. The 
equations include the von Mises-Hencky equation. 
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(S Q /S e ) 2 + (S B /S U ) Z =1 { 7 . 5 . 1 ) 

and the Gerber parabola equation 

(S /S ) + (S n /S u ) 2 = 1 (7.5.2) 

Where = alternating stress 
S = mean stress 
S e ~ endurance strength 
S u b ultimate tensile strength. 

In both of these eqxiations the information available from 
the Goodman fatigue diagrams presented in Chapters IV and 
V specify all of the quantities except the ultimate strength. 

A method was then sought which would give an estimate 
for the ultimate strength. It was assumed that the von 
Mises-Honchy equation was a valid mathematical model for the 
fatigue data. Graybill (9, p« 111) presents a method 
which can be used to calculate the least squares estimator of 
the ultimate strength assuming that the von Mises-Henchy 
equation adequately, models the fatigue data. The only 
other assumption which needs to be made is that the fatigue 
data to be used is not in the low cycle fatigue range, the 
low cycle fatigue range being below 10^ cycles. It was 
expected that the ultimate strength which would be predicted 
in this range would be quite large. This is quite a valid 
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assumption as low cycle life requires a completely different 
mathematical model than the von Hises-Honcky equation (See 
Chapter XI). 

The first step of the estimation process for the ultimate 
strength using the von Mises-Hcncky equation is to transform 


that equation into the following form 

1) (S Q /S e ) 2 + (S M /S u ) 2 =l (7.S.3) 

2 

2) Set y=S Q , and subtract (S a /S n ) 

from both sides yielding (y/S 0 )^-l - ( x / s u )^ (7,5,4) 

3) setting y'=S e ^-y^ (7.5.5) 

B " S e 2/S u 2 (7.5.6) 

4) Substitution yields y’-'B^x^) (7.5.7) 


The least squares estimate of B is given as 

»' = GO T GO H T [y] ^.s.7) 

where the brackets indicate vector quantities. Unfortunately, 
on Figure 7.5.1 indicates, there are only four values to be 
placed in the x vector. Ihesc values of the mean stress 
are derived from the y, or alternating stress values, in the 
follov?ing manner. Each of the mean stress values are related 
to the alternating stress value by the stress ratio r; where r 
equals 



(x 2 y 2 ) 




fXjO) 


Goodman Diagram Illustrating Mean of Strength 
Distributions Used in Least Squares Estimate 
of Ultimate Strength. (von-Mises Hencky Equation) 


>>>u^liiuiiUtUllMWj44UMl«M^tUUil 'HIM MUIUltU I :i< il Of U 1 1 IHJ‘l IIH ,1 
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x ~ s a / s m (7.5.8) 

It was possible to solve for the three values of mean stress 
by dividing the alternating stresses by r for each Goodman 
diagram investigated. 

Tims, the x vector was specified by three values of 
mean stress while the y vector was specified by the 
corresponding values of alternating stress. The value of 
B‘ is the least squares estimator of the X axis intercept 
which can then be related to the ultimate strength by 
Equation 7.5.6. 


Revrriting equation 7.5.7 in terms of the von Mises-Henchy 
transformed variables yields; 



It can readily be seen that x is zero which will cause 

0 

the vector equation directly above to be reduced to / 
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Hie first two vector multiplications t x ) f x I yield o scalar. 

The inverse of this scalar is simply the numerical inverse. 

Hie remaining vector multiplications are straight forward. 

The vector operations; after being transformed to algebraic 
relationships can be programmed (LSEFD) on the PDP-8 computer for 
the Goodman diagram data developed in Chapters IV and V of this report. 
A flow chart, variable definitions and computer listing of LSEFD pro- 
gram is given in Appendix E. 

The results of the least squares estimate is given in 
Tables 7.5.1 for Chapter IV Goodman diagrams and Table 7.5.2 
for data extracted from the Goodman diagrams of Chapter V. 

Cycles to failure data and the corresponding strength distribu- 
tions which this data would yield above the cycle life of 
200,000 cycles is unavailable at this time. It appeal's from 
Tables 7.5.1 and 7.5.2 that as cycle life increases the 
least squares estimate of the ultimate strength distribution 
decreases. Had the least squares estimate predicted consistently 
a value of 255,300 psi, for the measured ultimate strength of the 
Phase I grooved specimens, it could have heen concluded that 
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Table 7.5.1 Chapter XV Goodman Diagram Data Used for a 
Least Squares Estimate of the Ultimate 
Strength. 


Cycle Life 


Least Squares Estimate 
of Ultimate Strength 


20.000 - 40,000-' 

60.000 - 90,000 

90.000 - 200,000 


238, 3S5 
264,662 
222,661 


Table 7.5.2 Chapter V Goodman Diagram Data Used for a Least 
Squares Estimate of the Ultimate Strength. 

Cycle Life 

40.000 

90.000 

200,000 


Least Squares Estimate 
of Ultimate Strength 

262,533 

234,972 

216,190 


\ 
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the grooved geonetry specimen’s ultimate strength distribution 
was the proper distribution to use along the mean stress axis. 
The decrease in the least squares estimate of the ultimate 
strength as the cycle life increases senmingly shakes the pos- 
sibility that it is indeed the grooved specimens' distributions 
that should be used along the mean stress axis until the limit- 
ing case is investigated. This limiting- case being an infinite 
life Goodman diagram, kfaen the least squares estimate tech- 
nique is applied to the infinite life diagram (3, p. 75), pre- 
sented in Figure 4.4.1 which was developed by John Smith for 
the grooved specimens under discussions, the estimate of the 
ultimate strength is 222,661 psi. This is well above the value 
of the ungrooved ultimate strength of 178,000 psi, and is 87% 
of the ultimate strength of the grooved specimens ultimate of 
255,300 psi. 

It is also important to note that there is no data avail- 
able for strength distributions below a stress ratio of 0.44. 
Because of this the least squares estimate was bassed on four 
points with a stress ratio greater than or equal to 0.44. 

It is noted that even in the Goodman diagrams of Chapter V for 
cycle life values of 200,000 cycles, the largest invest5.gated 
alternating stress level at a stress ratio of 0.44 has only 
fallen 17% from the value at the stress ratio of infinity,, 
while the mean stress value is already 76 % of the ungrooved 
ultimate tensile strength at that particular point. 


The least 
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squares estimate would have been considerably more accurate 
if data below a stress ratio of 0.44- was available. Xt is in 
this region that the curve must transition to either the grooved, 
or the ungrooved ultimate tensile strength. Perhaps if this 
data were available the estimates would not have fallen off 
to the values below the grooved ultimate strength. The analytic 
results, which suffered from a lack of data below a stress 
ratio of 0.44, show that they reach a limiting case value much 
greater than that of the ungrooved ultimate strength, and con- 
firm the opinions which were solicited from noted authors on 
the subject of fatigue. It is the considered opinion of this 
investigator that the ultimate strength distribution of the 
grooved geometry test specimen should bo used for Goodman 
diagrams which represent the behavior of grooved fatigue test 
specimens. 


CHAPTER VIII 


EMPIRICAL MATH MODELING OP FINITE LIFE GOODMAN DIAGRAM 


8 . 1 Mathematical Models of the Goodman Diagram 

A mathematical model of the Goodman diagram relates this 
alternating stress/ S A' to the mean stress S m' by means of an 
algebraic equation. In most cases this equation relates a 
static strength parameter to the endurance strength, S e , of 
the material in question. The safe design region of a Goodman 
diagram is conventionally defined as the area bounded by the 
ordinate and abscissa axes and the line of the mathematical 
fatigue model, or equation, under consideration. There are 
several mathematical models of the Goodman diagram. These 
include the modified Goodman line, the Gerber parabola, the 
von Mises-Henchy ellipse, the Soderberg line, the Sines line, 
and the Longer modification to the modified Goodman line, The 
objective of this chapter is to present these mathematical 
models, however the presentation of such models would be in- 
complete if not accompanied by a discussion of the strengths 
and shortcomings of each model. Figure 8.1 compares the six 
mathematical models discussed in this Chapter, 

8 . 2 Modified Goodman Li no 

The most widely accepted theory of combined stresses is the 

modified Goodman line. The modified 
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Alternating Stress 
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Goodman line connects the endurance strength, S e , on the ordinate 
axis by a straight line to the ultimate tensile strength t s u , 
which is plotted on the abscissa axis, (10, p. 6). 

It is important to note that the endurance strength must be 
defined in relation to a given number of cycles, beyond which 
the material is assumed to have an infinite life, A common 
cycle life value for the endurance strength of ductile steels is 
10^ cycles. The equation of the modified Goodman line is: 


S </ S e * V S u = 1 


A common criticism of the Goodman line is that it tends 
to be conservative where the stress ratios, r - S G /S^, is 
well above one or in the range 1 < r < » . 

Even though it is conservative in this range, for lov/er stress 
ratios in the range of r = 1/10 and less, the modified Goodman 
diagram may predict safe combinations of alternating and mean 
stress when in actuality they could cause yielding (10, p. 7). 

8 • 3 Gerber Parabola - 

The Gerber Parabola was first proposed by Gerber in 1874, 

At this time Gerber was attempting to fit a curve to the resul.ts 
of Wohler's experiments with combined stresses. The parabola 
which Gerber proposed, known today as.. the Gerber parabola, joined 
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the ultimate tensile strength on the abscissa to the cn&urcmce 
strength on the ordinate axis. These tv/o physical properties 
uniquely describe the shape of the Gerber parabola. The Gerber 
parabola requires that a parabola be drawn so as to have a 
vertex on the vertical axis at the value of the endurance 
limit and to pass through the ultimate strength which is plotted 
on the horizontal axis (10, p, 8). The equation of the Gerber 
parabola is given by 

V S e + <V S u> 2 ° 1 

Critics of the modified Goodman bine have stated that it 
is too conservative. The Gerber parabola was proposed to com- 
pensate for the conservatism of the modified Goodman line . In 
addition, it has been found that in many cases the Gerber para- 
bola fits the experimental data far better, in the stress ratio 
ranges of 1 r ^ co , than the modified Goodman line. Unfortunate- 
ly/ the Gerber parabola does not give a proper representation 
of the fatigue data where stress ratios of one tenth and less 
arc encoitntercd. At these stress ratios the Gerber parabola 
permits an even greater amount of yielding than the modified 
Goodman line . 
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^ • 4 von Mtscs-Henchy Ellipse 

This equation, associated with the energy of distortion 
■theory discussed in the following chapter, has been proposed 
for the case of cor' ined stresses, Ibe von Mises-IIencky equa- 
tion or ellipse is given by 

(W z - <V S U >* = 1 

This equation forms an ellipse in the first quadrant of the 
cartesian plane. Although this equation was originally pro- 
posed for static loads, it is commonly used .as a model of com- 
bined stresses in fatigue (10, p. 1C), and (it is not theoretically 
valid above the yield point of the material). 

8 . 5 Sodorberg Line. 

In 1930 Scderberg proposed his theory in the United States 
which was to eliminate the problem discussed previously concern- 
ing yielding in the safe design region. The Soderbcrg line 
eliminates the problem of the yield point of the material being 
exceeded at any combination of stress. If yielding does occur, 
the dimensions of the specimen are changed. Obviously, this 
change is of a permanent nature and the performance of the 
material is affected. Even though failv.ro of the material nay 
be considered as fracture, the maximum allowable stress level 
becomes the yield strength. The actual Sodcrberg line tahes 
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these facts into account by specifying a straight lino between 
the endurance strength on the ordinate axis and +he yield 
strength Sy on the abscissa axis (11, 0,15). The equation of 
the Soderberg line is given by 

s „ - s . u - vy 

Criticism of the Soderberg line arises because it does 
indeed lie below the Goodman line. Because of this fact it 
will be even wore conservative than the Modified Goodman line 
for stress ratios in the range of 1-^r-^co, the Goodman line 
has been shown to be conservative in this region (11, p. 15, 16, 
p. IS). 

8 . 6 Sines bine 

The Sines line is an empirical relationship given by 

S = S - cS 
a e m 

where tho constant c must be determined for the material under 
investigation. Because it is an empirical relation it can 
through the appropriate value of c, be ad-justed to fit the 
data for a particular naterial. The Sines line accounts for 
maximum stresses up to the yield point. Consequently, the Sines 
line is defined only to the vertical lino where the nean stress 
is equal to the yield strength as shown in Figure 8.1 (10, p. 9). 
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8,7 Longer Modificat ion to the Modified Goodman Lino 

The Longer Modification to the modified Goodman Jino at- 
tempts to solve the problem of yielding for a different case 
than discussed in the Soderborg Line presentation. It will be 
recalled that the Sodcrberg lino attempted to solve the problem 
of yield at low stress ratios . The Longer modification attempts 
to. solve the problem of yield caused by high stress levels. At 
high stress levels the maximum val\ie of the alternating and 
mean stress. S + S M , which are encountered may exceed the 
yield strength (11, p. 16), As discussed in the Soderborg 
line presentation this yielding will have an adverse effect 
Upon the fatigue characteristics of the specimen, The Longer 
modification excludes the araa of the safe design region where 
the alternating stress plus the mean stress is greater than the 
yield strength. Hence, the safe design region becomes the area 
bounded by the ordinate and abscissa axes, the Goodman line and 
the region which satisfies the inequality S_ + S ^ S , 

. a m y 

Criticism is again directed at the Langer nodSfication 
because it is considered too conservative at stress ratios 
greater than one, The Langer modification docs eliminate the 
criticism of the modified Goodman line where the yield strength 
is exceeded and at the sane time does not demand the conserva- 
tism, of the Soderberg line (11, p. 16), 


CHAPTER IX 


THEORETICAL STRENGTH THEORIES 
9.1 Introduction to Strength Theories 

Ihe three principal stresses completely describe the 
stress state of any point in a structure. The need for a 
strength theory to describe the material behavior at a point 
arises when two or more of the principal stresses have a non- 
zero value, When one of the three principle stresses is non- 
zero the behavior of the mater5.al is described by the conven- 
tional tensile test. There are however/ even in this relative- 
ly simple state of stress , differences between the true state 
of stress and the "engineering” stress-strain properties which 
have been previously discussed in Chapter VII of this report. 
The objective of a "theory of strength" is to relate a complex 
state of stress/ i.e., v/hen two or more of the principal 
stresses are non-zero, to the uniaxial properties which are 
obtained in a tensile test (13/ p, 1). 

The elastic portion of the total strain is related to 
stress by Hooke's Lav;. In the case of combined stresses ini- 
tial yielding must be related to yielding in a tensile test by 
means of a flow theory. A flow theory relates the increments 
of plastic strain in each direction to the state of stress at 
the point under consideration. In the past there have been 
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several yield criteria, or* strength theories, proposed but 
because of later experiments in hydrostatic stresses, vrhi.ch 
have conflicted with those theories, they seen now to be only 
of historical importance (13, p. 1). However, there are tiro 
such theories which do not have this particular fault# These 
two are the maximum shear stress criterion, and the energy of 
distortion or von Mises-Henchy criterion. The following chapter 
Vfill discuss these two theorie.s, and will present three addi- 
tional modified strength theories proposed by Findley and 
Mathur (14). These discussions will be accompanied by compari- 
son of these theories to the fatigue problem of combined 
stresses of bending and torque. 

9,2 Energy of Distortion Theory 

The octahedral shear stress, the energy of distrotion, or 
the Von Hises-Henchy theory, as this theory is often referred 
to, predicts yielding to occur when the elastic energy of dis- 
tortion reaches a critical value. The energy of distortion is 
defined as the total energy minus the energy associated with a 
volumetric dilation. It can be shown that the energy of dis- 
torsion is proportional to the shear stress on the octahedral 
plane . The octahedral plane is the plane which nahes equal 
angles with the three principle directions (13, p. 3). 

Consider a cubic element of material acted xipon in the 
three principle directions by the stresses e^, Sp ( , and Sg 
where s^. For the unit cube, pictured in Figure 9.2.1, 








\ 


the work done in any principle direction is given by 


U -st/2 
n n n 


(9.2.1) 


where £ n are the three principle strains. Considering 


e - sj/E -/ts 2 /E “ 


(9.2.2a) 


c 2 = c 2 /E -^s^/E - AS 3 /E 


(9.2.2b) 


e 3 = s 3 /E - H s-^/E - M s^/E 


where At ~ Poissons ratio and E = Modulus of Elasticity. 


(9.2.2c) 


The total strain energy is 


2 2 2 


U = U 1 + U 2 -KJ 3 = 1/2E { Si +s 2 +c 3 - - 2JX(^ Z + c 2 s 3 + 838 ^- 


(9.2.3) 


Defining average stress as 


S 1 + s 2 + s 3 


(9.2.4) 


which is applied to each of the principle directions of the 

unit cube, the remaining stresses s, - s , s_ - s and 

x avg' 2 avg 

s 3 - s , shown in Figure 9.2, je will only produce angular 

distortion. If s is substituted for s. . s„ and s 0 in Equa- 
avg 12 3 A 

tion 9.2.3 the amount of strain energy which produces only change 


in volume it 


! 




132 


U •» ( 3s ova - 2 -H(3 a avg> > ={3^ vg /2E X1-2M) (9.2.5) 


2K 


avg 


“ { (s, + s, + So ) } 


3 


If the expression for is substituted in Equation (9.2.5) 

it becomes 


2 2 2 

= {1-2M/6E }(s^ +S 2 +Sg + 2s2Sg + 2f3 3 s i) 

(9.2.6) 

The energy of distortion is then equal to the total energy, 
given by Equation 9.2.3 minus the energy of the volume change 
givc-n by Equation 9.2.6. The energy of distortion is thus 
given by 

2 

U d = U - U v =(l+^'3E){(8 1 .sp j ) 2 +(6 2 -B 3 ) 2 +(s3-s ;L ) 2 } 

V/hen a state of pure shear exists, the shearing stress 
at point, ‘‘t , is equal in magnitude to each of the principle 
Stresses at the same point. If = 't and S 2 ~~ s y the energy 
of distortion becomes 

2 

U d = (lW<SE) {(s-(-s)) + (s) 2 + (-s) 2 } *{1+ ^s 2 /S} S 2 

V/hen s ~X , then ' 

G 2 

=(l+^VE)r c . 


which is the energy of distortion in p torsional test specimen. 





133 

If the sane unit cube is subjected to a normal stress, s^, 
in one direction only, the other two principal stresses being 
zero as in a tensile test spec.inen subjected to an axial load, 
the energy' of distortion, U^, becones 

2 y 2 

U d = {1+ /V6 E + Sj ') = CL+*V3E )b 1 

Equating the energy of distortion for the case of pure 
shear to that of the unaxial tension condition, when yielding 
first occurs ; 

(H- /-Vk) = (l+"-/3E) n Q 2 

t 0 = S e /V^ 

The conclusion which can be reached is that yielding and ' j ' - 

eventually a ductile fracture starts when the energy of dis- 
tortion reaches a critical value. The maximum shear stress, 

Z c , a point when yielding starts is 1/ 3 tines the 

maximum tensile stress, s e , at the sane point. Thus the 
von Mises-Henchy ellipse for combined bending and shear stress 
is given by | • 

. (V s .) 2 ♦ <V S „>* = 1 i' 

(W* ♦ (0/VTH) 2 - 1 j , 

I I 

where in the above equation is the constant shear stress 
recorded in each test specimen studied. 
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9.3 Maximum Shear Stress Theory 

Tresca proposed in his maximum shear stress theory that 
yielding occurs when the maximum shear stress reaches a criti- 
cal value (13, p. 2). Yielding begins when the maximum shear- 
stress equals the shear. stress corresponding to the yield strength 
in the simple tension test. According to Figure 9.3.1 yielding 
occurs when X - S /2 where S„. is the yield strength of the 

iiiCvv y y 

material. For a triaxial stress state three maximum shear 
stress may be found and are given by 

X =(s i - s^/2 X = (s 2 - S 3}/ 2 X = (s JL - s 3 )/2 

Yielding vrill begin when the largest of these shearing stresses 
becomes equal to one-ha 1 f the tensile yield strength of a simple 
tension test specimen (15, p. 152). 

In essence the theory predicts that the shearing yield 
strength is equal to one -half of the tensile yield strength. 

The advantages of the theory is that it is easy to use, is 
useful for ductile materials and is conservative in describing 
the behavior of brittle materials (15, p. 152). 

9.4 Comparison of the Maximum Shear Stress and Energ y 
of Distortion Theories 

The maximum shear stress and energy of distortion theories 
can be conveniently represented in Figure 9.4.1 in the two- 
dimensional principal stress space. 
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Fig. 9.41 Representation of Energy of Distortion and 
Maximum Shear Stress Theories in the Plane 
Perpendicular to the Unit Vector (l/iTT*, l/VT", l/fT' ) 
Cl 3 1 p, 81 
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Experimental evidence indicates that the Von Mises-Henchy 
or energy of distortion, criterion is more accurate than the 
'fresco theory for predicting the yield strength of most material 
under biaxial stress. In addition strain hardening and creep 
behavior correlate much better for most materials using the 
on Misee-Henchy theory. The difference between the two theories 
is small, the maximum difference in the two theories in any 
state of stress being about 16 percent. In view of this fact the Tresca 
theory, v:hich is more conservative is considered satisfactory 
even though the von Mises-Menchy theory is more accurate 
(13, p. 6). 

It has been stated that stress -strain properties are best 
correlated using the von Kisos-Henclcy theory. This, however, is 
not grounds enough to state that fatigue failures are best 
described by the same criterion. The energy, of distortion 
has no directional properties and is always considered a posi- 
tive quantity. This causes some serious deficiencies as a 
means of predicting fatigue failure. In fatigue experiments 
it has been found that the maximum shear stress theory will 
correlate results as well as the von Mises-Henchy criterion. 

'ilia problem of determining which of the' theories is best is 
difficult because of the natural scatter of fatigue data spanning 
the various criteria. 

Referring again to directional properties of the two 
criteria it is noted that the maximum shear stress changes sign 
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when the stresses ore reversed, os during a rotciting bectra 
fatigue test, whereas, it has been previously noted that the 
energy of distortion is always a positive quantity. In re- 
versed bending tests these directional properties are of crit- 
ical importance. As the stresses at a point are reversed the 
energy of distortion goes from a positive quantity to zero and 
then bach to a positive quantity. This means that the distor- 
tion energy theory is not reflecting a reversal in loading. 
Reversed loading cases are very harmful to the fatigue of any 
structure. It is important to recognize load reversals as 
adding to the range of stress and strain. It would be possible 
to devise methods to account for load reversals in the uniaxial 
case, however, in the case of combined stress the problem would 
be quite complex (13, p. C). 

The von liiees-Henchy theory is not free of the directional 
property problem in fatigue. As an example one can consider 
the special case where the magnitudes of the principal stresses 
are constant but the directions arc changing with time. The 
energy of distortion remains constant indicating that there is 
no fatigue loading. However, this is not the case as the shear 
stress theory must also be modified for the case where the direc- 
tions of the principal stresses vary with time. However, in 
this case one would merely use the history of maximum shear 
stress on a fixed plane of the material to predict the fatigue 


life (13, p. 9). 


Finclley and Kalhur (14) have developed several theories of 
fatigue failure under combined stresses. The following section 
will present the three theories which were developed by then. 

They state that the classical theories for the initiation of 
yielding under combined stresses are conflicting. The investi- 
gation of Gough and Follard (14, p. 2) concltided that the 
classical theories which have been proposed are inaccurate 

since the ratio of fatigue strength in bending, b “ S q / to 

°e 

that in torsion fatigue strength, t = S^, is not the sane for 
all n eta Is as required by the classical theories of yielding. 

Findley and Hathur propose that the. influence of the prop- 
erty of anisotropy, that of having different properties in 
different directions, and the state of combined stress are the 
cause of discrepancies between proposed theories of failure 
o.nd results obtained from combined stress fatigue tests (4, p. 3) 

The ratio of b/t varies over a considerable range for all 
engineering materials. Considering only metals this ratio varies 
from a value of 0.9 to 2.6. If, however, the metals are grouped, 
within each group the ratio of b/t ha3 a much smaller range. 

Cast irons and its alloys have a b/t ratio which varies from 1.3 
to 2.5 with a majority of these values lying between 1.5 and 
2.0, This latter range is considered to be the ductile range. 

The b/t ratio value cf notched steels was found to be consider- 
ably less than that of unnotched steels with a majority falling 
between 1.0 end 1,5. 
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It is felt that anisotropy of the material is a contribu- 
ting caxjse to the variation of the b/t ratio. It has been 
suggested that fatigue failure is caused primarily by an alter- 
nating shearing stress producing repeated slip, however, the 
resistance of a material to this fatigue mechanism may be 
influenced by the magnitude and sign of the normal stress oc- 
curing on plar.es of maximum shearing stress. This effect 
may vary with the material. Differences in the effect which 
the normal stress may have in fatigue in a given material will 
cause differences in the b/t ratio (14, p. 6), 

9.5.1 Correction Factors 

In ductile metals the. cyclic principle shear stress 
is the quantity most closely associated with tire fatigue damage . 
'fire principle shear stress theory predicts that the bending 
strength should be twice the shear strength (14, p. 8). How- 
ever, because of the effect of anisotropy and combined stress, 
Findley and Hathur suggest that the principle shear stress 
theory be modified to the form 


max 


= b/ZK s t or K - b/2t 


(9. 5. 1.1) 


where K is the correction factor for anisotropy and the combined 
stress condition. The expression for modified principal shear 



(9. 5. 1.2) 





Where end represent respectively the amplitudes of alter- 
nating bending stress and alternating torsional stress com- 
ponents of combined stress (14, p, 9), 

It is important to note that v/hen correction factors are 
applied to the octahedral shear stress, principal shear strain, 
energy of distortion, total energy of distortion, total energy 
of deformation and magnitude of state of stress vector 
theories the result is the same in all cases, that of Equation 
9. 5. 1.2. 

If correction factors are applied to the principle stress 
theory the governing equation for combined bending and torsional 
stress becomes 

V S S E s n /s u = 1 (9. 5. 1.3) 


When the principle strain theory is modified the fatigue 
strength ratio b/t predicted is 


/S - b/t = 
„ u 


1 + p. 


(9. 5. 1.4) 


where JJ. is Poisson's ratio. 

The expression for the modified principal strain theory 
as reported previously (17) is 


~ 1+ M /2r, +f2b/t(l+/<A.|r + (l- A</l+/c) s| 

For pure bending 

% “ 0, S = b, thus i i " b/E 


(9. 5. 1.5) 


142 


Substituting for e in equation 9 . 5 . 1 . 5, and 

JL 

simplifying the expression, the modified principal strain 
theory becomes 

JUL s 2 /b 2 + (1 -/*- ) s/b +bV “ 1 (9. 5. 1.6) 

The fatigue strength ratio b/t predicted by the principal 
strain theory (17) is 

b/t = 1 + JLL (9. 5. 1.7) 

Substituting equation 9. 5. 1.7 in equation 9. 5. 1.6 

(b/t - 1) s 2 /b 2 + (2 - b/t) s/b t 'i/ 1 2 a 1 

(9. 5. 1.8) 

This equation models the combined stress state of bending and 
torsion for the modified principle strain theory. 

A design expression has been px'oposed by Findley and Hathur 
(14) to model the fatigue failure of notched ductile metals and 
irons. These materials have a behavior which is intermediate 
between the perfectly brittle irons which have little or no 
slip and ductile metals which have considerable slip. The 
stress system which is associated with these materials may 
change from the principle stress theory, b/t = 1, for materials 
such as the brittle irons to the principle shear stress theory 
where b/t - 2.0 as in the ductile metals (14, p, 11). 


14 


The design expression which is suggested to model these 
metals under combined bending and shear is given by 


( VV 


b/t 2 

+ (s /s u ) = >• 


m 


(9. 5. 1.9) 


It is noted that the exponent b/t varies with, the class of material . 
l/hen b/t = 1 Ecpacition 9. 5, 1.9 is that of the modified principal 
stress theory. V/hen b/t “ 2, vrhich is the case in ductile 
materials the equation reduces to that of the modified shear 
stress theory (1 4, p. 11). 

9,5.2 Comparison to Faticrue Data 

Findley and Hahur compare their modified stress 
theories to actual combined bending and torsional fatigue test 
results. It is found that the modified sheaf stress theory. 
Equation 9. 5. 2.1 served as a good model for ductile metals with 
a b/t ratio ranging from 1.46 to 2,0. A comparison of this 
equation to actual fatigue tests is given in Figure 9.5.1 
(14, p. 15). 

The modified principles stress theory is compared to actual 
fatigue data of iron and iron alloys in Figure 9.5,2 and to thet 
of notched ductile steels in the brittle range, b/t = 1.3 in 
Figure 9.5.3. The modified principle strain theory predicts 
strengths which are higher than the actual data as can be seen 

















Figure .9. 5. .3 Comparison of Design E 
Under Bending and Torsional Stre 
Ratio Greater Than 1.3. (14, • 








Figure 9.5.4 Position of Modified Strength Theories as Specified Bv 
b/t Ratio. (14, p'. 25) ‘ 

Curve A: Principle Stress 'Theory 
Curve B: Principle Strain Theory 
Curve C: Total Strain Energy Theory 

Curve D: Distortion Energy Theory 

Curve E: Principle Shear Stress Theory 
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in Figure 9,5.?. (14, p. 15). 

It is found that the design expression. Equation 9. 5. 1.9 
is in very good agreement with data ox notched ductile steels, 
b/t * 1.3, and in good agreement with both brittle and ductile 
range of metals as shown in Figure 9.5,3 (9, p. 15), 

Assigning the metals to regions by b/t ratio values, as 
discussed previously, yields the following regions: 

Region I (brittle) 

Region II (intermediate) 

Region III (ductile) 

The position of each modified strength theory within the regions 
described in Figure 9,5.4 is specified by the value predicted 
by the b/t ratio of each strength theory. Table 9.5.1 presents 
the value predicted by each strength theory for the ratio of 
b/t. The modified strength theories may be assigned to each 
of the regions described previously in the following descending 
order (9, p. 17). 

Region I Modified principal stress 

Design expression 
Modified principle strain 

Region II Design expression (notched ductile steels) 
(equation 9. 5. 1.9) 

Modified principal strain 
Modified principle stress 
Modified principle shear stress 

Region III Modified principle shear stress 

Design expression (equation 9. 5, 1,9) 


Table 9.5.1 Values of the Ratio b/t Predicted by each Strength 
Theory (14, p. 17). 


Principle stress theory 

Principle strain theory 

Total strain energy theory 

Distortion energy and 

octahedi'al shear stress theory 

Principle shear stress theory. 


b/t *= 1 
. b/t « 1 + P- 
b/t + 2 M- 

b/t - Vf 

b/t =2.0 


CHAPTER X 


RECOMMENDED EMPIRICAL MATHEMATICAL MODELS OP THE 
FINITE LIFE GOODMAN DIAGRAM 

In viewing each of the five finite life Goodman diagrams 
developed in Chapter V it can bo seen that tho equation of 
the mean line is generally of a quadratic nature. It is 
of interest to investigate if the mean line of the finite life 
Goodman diagrams can be described by one or more of the empirical 
mathematical models of the Goodman diagram discussed in Chapter 
VIII. 

The Von Mises-IIenclcy ellipse was compared to the experi- 
mental finite life Goodman diagrams presented, in Figure 5,3,1 
throxigh 5.3,5, The ellipse is given by 

<S Q /S e )2 ♦ (S„/S u ) 2 = i (10.1) 

The endurance strength and ultimate strength mean values speci- 
fy the mean line of the von Kises-Hencky ellipse. The enduro.nc 
strengths for each of the five von Mises-Hencky ellipses were 
taken from Table 5.3,1, These ellipses, superimposed upon the 
original experimental mean line, appear as overlays with the 
experimental finite life Goodman diagram in Chapter V as 
ove slays on Figures 5,3,1 through 5.3,5 As in Chapter V the 
ultimate strength distribution is taken to be that of the 
grooved specimen, A PDP-8 program, described in Appendix E, was 


150 


151 


used to calculate the mean stress points as specified by the 
ellipse for several theoretical alternating stress levels. 

Those values are listed in Table 10,1 for each of the five 
finite life Goodman diagrams. 

The finite life Goodman diagrams for 40,000, 90,000, and 
200,000 cycles are in very close agreement with the von J'ises- 
HencJcy ellipse. The lower cycle life diagrams of 3,500 and 9,000 
cycles are in very poor agreement with the von Mi.ses~]Jencky 
ellipse, a fact which will be investigated later in this 
Chapter, 

Although the higher cycle life fatigue diagrams seemed to 
be closely approximated by the von Kises-Henchy ellipse, it was 
decided to determine the exponent, a, where in general the 
quadratic equation of interest is given by 

< S a /S e>° + (W* “ 1 (10-2) 

This was accomplished in the following manner. Equation 10.2 
can be transformed to 


a 7. 

1 + x “ 


(10.3) 


y ~ S /s 

J a e 


X ~ S /S 
m u 


where 


0 

0 


I 



Tabic 10.1 Alternating anti Mean Stress Values Predicted by the von 
Mises-Hencky Ellipse. 


S - psi S - psi 

a . in 


N=3,500 cycles; 

134,000 

43,118 


128,000 

86,038 

S -135,953 psi; 

120,000 

119,995 


100,000 

172,959 

S u =255,300 psi 

t>U, GOG 

2^y,UW2 

20,000 

252,522 


N-9,000 cycles; 

S =115,926 
e 

112,000 

168,000 

100,000 

80,000 

60,000 | 
20,000 

65,878 

92,779 

129.145 

184.145 
218,445 
251,472 

N=40,000 cycles; 

88,000 

61,751 



84 ,G0C 

96,000 


S =90,693 

80,000 

120,264 


c 

60,000 

191,445 



40,000 

229,128 



20,000 

249,015 


N=90,000 cycles; 

68,000 

130,204 



64,000 

149,866 


S =79,054 

60,000 

166,231 


e 

40,000 

220,207 



20,000 

246,995 


N=200,000 cycles; 

68,000 

63,027 



64,000 

104,696 


S =70,172 

60,000 

132,388 



40,000 

209,761 


. 

20,000 

244,711 



UrnWri***!* 


153 


£ 

Solving for y yields 

y a = 1 - x 2 (10.4) 

If the natural logarithms of both sides are taken, the result is 

a In y = In (1-x 2 ) (10.5) 

In y = 1/a In (1-x 2 ) (10.5) 

If this equation is plotted on In - In graph paper the result will be a 

straight line. A least squares PDP-8 program SBFB, is available which 

will give the slope of the best fit equation through the data as well 

as the correlation coefficient. The slope of this line is equal to the 

value of the unverse of the exponent a in Equation 10.2. A description 

of this program is given in Appendix F.. The values derived by this 

method for the exponent a and p are given in Table 10.2 

The finite life Goodman diagrams of 3,500 and 9,009 cycles as 

mentioned previously of not seem to be modeled by a quadratic equation. 

The various empirical models were reviewed and it was deterimed that 

the Sines Line had the greatest potential of modeling the finite life 
4 

range below 10 cycles. The Sines Line is given by 

m 


S - S - cS 
a e i 


( 10 . 6 ) 




Table 1° . 2 


Values of the Exponent a and Correlation Coefficient. 

Cp) 


Cycle Life-N 

Coeff. a 

Correlation Coeff.(p) 

1*0,000 

2.3821 

• 9137 

i 

90,000 

2.3042 

.8948 



200 >000 


1.9107 


.8915 
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The coefficient c for the finite life Goodman diagrams of 3 / 500 
and 9,000 cycles was found, by trial and error, to be 0,03, 

The Sines Line specified by c - 0,00 v/as placed on overlays of 
Figures 5,3.1 and 5.3.2. The Sines Line is valid only to the yield 
strength of the material. Because the yield strength of the 
grooved geometry specimen was not available, a dashed line v/as 
drawn beyond the yield strength cf the ungrooved specimen. 

This is meant to signify that the yield strength of the grooved 
specimen is some value greater than that of the ungroovod. speci- 


men. 

In conclusion, it can be said that the von Mises-Henchy 
ellipse quite closely approximates the mean line of the finite 
life Goodman diagram, where the cycle life value is above 10 r 
cycles. If the Von Hises-Kencky ellipse is modified to 


( s Ai ! V ) 2 - (V s .. i^s > 2 = 1 

e u 


(10.7) 


the Goodman surface is adequately described. This surface is 

shown in Fipure 5.3.5. For life values below 1C 4 

cycles the Sines Line adequately describes the mean line of 

the finite life Goodman diagrams. If the Sines Line is modified 

as shown in overlays of Figures 5,3. land 5.3,2 to the form 


S „ = < s e ± 3 °s ) " cS m / 


a 


( 10 . 8 ) 


then the Goodman surface is completely modeled. 
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The comparison of the von Mises-Hencky ellipse with the experi- 
mental finite life diagrams, and the discussion of Section 9.3 leave 
little doubt that the von Mises-Hencky strength theory closely pre- 
dicts the fatigue failure in steel. The discussion of Section 9.3 
gives adequate documentation that the von Mises-Hencky strength 
theory is widely considered as a strength theory for combined stress 
conditions. Although the von Mises-Hencky strength theory is not 
free of all criticisms when considered as the criterion governing 
fatigue failures the very close behavior of the larger cycle finite 
life Goodman diagrams indicate that this criterion can be continued 
to be accepted as governing the fatigue data generated under National 
Aeronautics and Space Administration Contract NGR 03-002-044 at The 
University of Arizona. The material which is used in this research 
program is SAE 4340 steel, Rockwell C 35/40. 


CHAPTER XX 


wo rscomiieiided methods of reducing the quantity 
EXPERIMKI iTAL DATA NEEDED FOR A FATIGUE 
DATA ACQUISITION PROGRAM 


OF 


The problem of reducing the number of data points in a 
fatigue test program is quite complex. The nature of the 
fatigue nechanism requires that distributions be developed to 
describe the cycles to failure and strength parameters of . 
the test specimen. To date the complex fatigue test program 
at The University of Arizona has performed fatigue tests on well 
over 650 specimens. This has been expensive but the experience 
gained by this pioneering program has laid the foundation to 
develop theoretical concepts vrhich may possibly reduce the need 
for expensive fatigue acquisition programs in the future , An 
examination of the Goodman diagrams and surfaces developed 
in this paper and the application of the theories proposed by 
Y7« N. Findley (14, pp. 2d) seen to offer great hope for reducing 
the number of test specimens required in the fatigue data acquisition 


program. 

In reviewing the finite life Goodman diagrams developed 
here, it becomes evident that the distributions placed along 
the mean and alternating stress axes are of critical importance 
in determining the shape ctnd location of the finite life 
Goodman diagram and surface. Currently the ultimate tensile 
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strength distribution which is placed along the mean stress 
axis is developed from the results of tensile tests. These 
tests include the pulling of thirty five grooved test specimens 
to determine the ultimate strength distribution; this method is 
satisfactory and should be continued. 

The distribution which is placed along the alternating 
stress axis can be determined as it is currently being done 
from the stress cycles to failure data at a stress ratio of 
infinity. This requires the generation of cycles to failure 
distributions at five alternating stress levels. Currently, 
this would require thirty five test specimens in each of the 
five stress levels tested. The methods presented in Chapters 
III and V could then he used to place the derived strength 
distributions on the finite life Goodman diagrams* Considering 
the conclusions reached in Chapter X the development of the 
finite life Goodman diagrams would be possible with two 
hundred test specimens, To review, the endurance strength foz* 
each cycle life would be specified by the transformation of the 
cycles to failure distributions at a stress ratio of infinity 
to the vertical strength distributions on the Goodman diagram. 
This procedure is outlined in Chapter V. The ultimate strength 
distribution is specified by the tensile test of the thirty-five 
grooved specimens. The entire finite life Goodman 

diagram could be modeled, above lc/ T cycles, by the Von IJises- 
Henchy ellipse as specified by/ 
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(S /S o + 3 c ) + (S /S '+ 3cb ) =1 (11.1 ) 

The necessary data would amount to a maximum of two hundred 
test specimens which would include a staircase analysis of the 
endurance strength at 10^ cycles, at a stress ratio of infinity. 

The required number of test specimens can bo reduced 
significantly if the cycle life values for the finite life 
Goodman diagrams are determined prior to the start of the test 
program* If for* instance, three finite life Goodman diagrams 
are required then the required endurance strengths can be 
determined by the staircase method for each cycle life. 

This method would reduce the number of specimens to one 
hundred twenty five to one hundred fifty depending upon the 
number of required diagrams. Each staircase method should have 
ci sample size of thirty five while the ultimate strength dis- 
tribution of the grooved specimen should also be specified by 
thirty five specimens . 


By drawing a desired stress ratio lines on the finite life 
Goodman diagram, as specified by the above equation, ,\i would 
be possible to develop the vertical strength distributions 
vhich would be placed on the statistical S-H diagram by follow- 
ing the reverse of the procedure developed in Chapter V. That 
is by drawing the stress ratio line on the finite life Goodman 
diagrams, the alternating stress level and standard deviation, 
in actuality the three sigma limits, of the rotated vertical 
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strength distributions would be obtained. It vwuld, by the 
reverse of the method discussed in Chapter V, be possible to 
create the vertical strength distributions which would be 
placed on the statistical C-N diagram. This, at least suggests 
that it may be possible to continue one step further back to de- 
termine the cycles to failure distributions. This 

method eliminates the need of generating cycles to failure dis- 
tributions at the intermediate stress ratios between zero and 
infinity. Currently such investigations require slightly more 
than four hundred specimens. This number considers thirty five 
specimens to a cycles to failure distribution. 

In an attempt to reduce the need for experimental data 
beyond that discussed above it appears that the theories proposed 
by Findley (14) in Chapter X offer an alternative to the expensive 
fatigue data acquisition program. Considering that the mean 
line of the finite life Goodman diagram follows the equation 

(S a /S e ) a (S/Su) 2 » 1 (11.2 ) 


The exponent a, which has previously been discussed in Chapter 
X also eppeared in the design expression of Equation 9. 5.1,9 
proposed by Findley for notched ductile steels. 


, . .b/t 2 . 2 

(s /s ) + s /s 

a c 151 vi 
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Although the equations 9. 5. 1.9 and 12.2 are not exactly 
identical differing in the denominators of the terms on the 
left hand side of the equations, they do attempt to relate 
quantities of bending and torsional stress , • Comparing the two 
equations, it can be seen that b is actually S e and t is S u . 

Hie exponent a is the ratio of bending strength to torsional 
strength, and could be determined experimentally. It would have 
to be determined if this ratio should be found from static 
tests or dynamic fatigue tests. In either event, facilities 
at Hie University of Arizona, including the NASA complex fatigue 
machines, v/ould.be adequate to determine this ratio, Findley 
Coleman, and Hanley site the value of b/t ~ 1.78 for an SAE 
4340 steel, Rockwell C 35 (13, p; 153); however, the dimensions 
of the test specimen are not the same as that undergoing 
tests at The University of Arizona. In addition the specimen 
which Findley, et. al., used in their studies is of the un- 
grooved geometry. It seems that the exponent a is in actuality 
Findley's exponent b/t. It would be possible to experimentally 
determine this value with no more thm the number of specimens 
required by the static ultimate strength tests, or thirty five 
specimens. 

Once the mean line is specified the standard deviation 
along this line could quite possibly be approximated by dividing 
the standard deviation of the endurance strength which is placed 
along the alternating stress axis by sin where © is specified 


162 


by the ratio of alternating stress to mean stress (3, p. 73), 

Hence, the standard deviation along the stress ratio axis, s , 

5 r 

is given by 

O - e O" /sin© (11.3) 

S S 
r e 

This then would give an approximation of the standard deviation 
of the Goodman surface, 

Iho two methods presented in this chapter, to reduce to 
a minimum the required number of test specimens, are not meant 
to be finalized proposals for a fatigue test program. They do 
seem to this investigator to be valid means, of reducing the need 
of large quantities of experimental data. The later method 
discussed would require only the number of specimens needed 
to experimentally determine the ratio of b/t, quite possibly 
no more than 35 - 40 specimens. The method formerly discussed 
v;ould require approximately 200 test specimens. Further in- 
vestigation of these two methods seems prudent,. 





CHAPTER XII 
OVERALL CONCLUSIONS 

Methodologies for developing finite life Goodman diagrams 
and surfaces have beer presented in this report. The finite 
life Goodman diagram presents allowable combinations of alter- 
nating and moan stress for the combined stress condition of 
alternating bending and constant torsional stresses for specific 
periods of design life. The actual Goodman surface, which is 
developed using two to five strength distributions at specified 
stress ratios, can be used to construct strength distributions 
at any desired stress ratio. The strength distributions are 
distributed normally. The technique of constructing a strength 
distribution at any specified stress ratio is initiated by 
construction of the desired stress ratio line on the Goodman 
diagram. The intercept of this line with the Goodman surface 
specifies the mean and the stendard deviation, in terms of 
the three signa limits, of the strength distribution at that 
stress ratio. The alternating stress level of the mean of 
the strength distribution can be read directly from the finite 
life Goodman diagram. This procedure is particularly valuable 
where the strength distribution at a specific stress ratio 
is required by the interference technique used in probabilistic 
design, 

16 * 
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In conclusion the finite life Goodman diagrams and surfaces 
developed in this paper indicate that as the cycle life decreases 
the allowable combinations of bending and shear stress magnitudes 
increase. The strength distribution to bo placed on the mean 
stress axis of the finite life Goodman diagram is concluded to 
be that of the same geometry test specimen as underwent fatigue 
tests. The grooved geometry test specimen is concluded to be 
the test specimen which specifies the ultimate strength distri- 
bution which is placed on the mean stress axis of the finite 
life Goodman diagram. The von Hises»UencJcy strength theory and 
ellipse have been shown to adequately model the behavior of 
SAE 4340 steel , above the cycle life of 10^ cycles, under 
combined alternating bending and constant torsional stresses. 


• CHAPTER XIII 


OVERALL RECOI D ISHDATIO! IS 

The follov/ing recommendations are offered by this inves- 
tigator: 

1* As additional cycles to failure data becomes available 
from Phase II of the national Aeronautics and Space 
Administration Grant No. 03-002-044 at The University 
of Arizona the methods $>roposed in Chapter XII to 
reduce the required number of test specimens should 
be further investigated. The recovery of cycles to 
failure data from the finite life Goodman diagram, 
should be investigated as proposed in Chapter XII , 
using both Phase I and Phase II data. 

2. A comprehensive literature search should be undertaken 
to obtain the complete set of papers authored by 
Professor W. U, Findley of Brown University. These 
works would be a valuable aid to the research program 
being conducted at The University of Arizona for the 
National Aeronautics and Space Administration. 

3. Additional fatigue and static strength data, beyond 
that supplied by Phase II of the National Aeronautics 
and Space Administration research effort, should be 
acquired through computer search facilities. These 
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facilities utilise high speed computers to selectively 
retrieve and display requested information , This 
would be particularly helpful in determining the 
best theoretical distribution to be. assigned to the 
static strength parameters* The additional data sup- 
plied by these search facilities would compliment 
data acquired through the experimental test programs. 
These search facilities include the Mechanical 
Properties Data Center in Traverse City, Michigan 
and the Defense Metals Information in Columbus, Ohio. 


APPENDIX A 


FORTRAN Computer Program To Reduce Cycles To Failure Data 


A-l. Flow Chart 

A-2. Definition of Variables 


A-3. Program Listing 
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MAIN PROGRAM 



Figure A-l. Flow Chart for Failure Data Distribution Determination Computer 
Program. 







Call subroutine for 
Chi-Square distribution 

flt ckifit) 


Print CHISQ , MCL 


GO T O 5 


END 


SUBROUTINE PROGRAM 


Subroutine KIFIT(N,) 
XN, XMEAN, XSIG 
CHISQ, MCL 

Calculate MCL, XMAX, 
XMIN, WIDTH, X(I) 
XMID(I) , OFREQ(I), 
EFROQ (I) 


C Print XN (U) , 

EFREQ(I) , OFREQ(I) 




Figure A-l (Cont'd) . 




Print reduced histogram 
EFREQ(I), OFREQ(I) , 

XCL 


U ... . "S 

” 

Calculate (x 2 ) CHISQ 

, 

... 

return ■ 1 

^ 

* 

r 

El 



Figure A-l (Cont'd). 



CONF s= TWO SIDED CONFIDENCE LEVEL 

T = STUDENT 1 FOR DESIRED CONFIDENCE M AND DEGREE 
OF FREEDOM N 

AREA = A3SISSA VALUE OF NORMAL DISTRIBUTION FOR 
SPECIFIED CONFIDENCE 

CHIU. CHI L = CHI SQUARE FOR CONFIDENCE M AND DEGREE 
OF FREEDOM ti 

S = STRESS LEVEL OF TESTED SPECIMENS 
N - NUMBER OF SPECIMENS 
XN ( 1 J s CYCLES TO FAILURE OF SPECIMENS 
XNLOG(l) = LOG CYCLES TO FAILURE OF SPECIMENS 
XNSUM = SUMMATION OF XNLOG(J) 

XKLOG = MEAN LOG CYCLES AT STRENGTH. S 

XSIGL = STANDARD DEVIATION IN LOG CYCLES AT STRENGTH 

XSIG = SORTF WITH K-l DEGREES OF FREEDOM 

XK2.XM3.XMA = 2ND, 3RD, 4TH MOMENTS OF LOG CYCLES 

XK3= COEFFICIENT OF SKEWNESS 

XK4 = COEFFICIENT OF KURTOSIS 

XMl.XMU = MEAN LOWER AND UPPER CONFIDENCE LIMITS OF 
MEAN LOG VALUE 

XSTDL.XSTDU - STANDARD DEVIATION LOWER AND UPPER 
CONFIDENCE LIMITS OF LOG VALUE 
CHI SO - VALUE OF Chi SQUARE F 1 T OF DATA HISTOGRAM TO 
NORMAL DISTRIBUTION DATA 

MCL = NUMBER OF HISTOGRAM CLASSES FOR CHI SQUARE FIT 
ADDITIONAL SUBROUTINE VARIABLES 
N = NUMBER OF DATA POINTS 
XKti) = DATA POINT VALUES 

XMEAN » XSIG --- MEAN AND STANDARD DEVIATION OF DATA 
VALUES 

xmax - maximum failure value 
XX in = minimum FAILURE value 

•WIDTH = CLASS INTERVAL WIDTH 
X f J ) ' CLASS BEGINNING POINT 
XMID(I) = CLASS MID-POINT 

OFREQ(I) = OBSERVED FREQUENCY OF OCCURRENCE IN CLASS 
INTERVAL I 

EFREQU) = EXPECTED FREQUENCY OF OCCURRENCE 


Figure A-2 - Definition of Variables for Failure Data Distribution 
Determination Computer Program. 
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*** SV-iTH.^b' 

* COMPILE FORT RA U» EXECUTE FORTRAN 

SLJSRO'JT I i\L K. 1 F ] T C ! i , X-\, X.MtiAN . XSIG, CM 1 SQ.KCL ) 

C CHI SOl'ARLD F i 1 OF DATA POINTS TO A NORMAL 0 1 ST 
46 FORMAT (lXiJlKLxPO'ficNT = .FS.O.l&K CLASS MIDPOINT 
1 FS.O) 

4 7 FORMAT (J 6 M OAT,\ POINT VALUES/ ( IX . FU . 3 .9F12. 3 ) ) 
t,i: FORMAT (11X.I0K EXPECTED » 1*0 X » 1 OH OBSERVED ,10X» 

1 7H CLASS/ J4X, 20HCL4SS FREQ. 1CX. 1CKCLASS FREQ * 
■2 ?>:, 3 ZHCEblNNiNG PT / (1 X , f 20 . S » F 20 . 1 » F 20 . 8 J J 

DIMENSION X.\(I 2 C). X(15)» Xv-lDdA). OF RFC < 1 5 ) , 

1 EFREcnSI 

Af; s N 



CLASS 

= -2 » C + 3 • 3 

*LOGF 

(AN) 


KCL * 

CLASS 




X:\AX = 

XN ( 1 ) 




X!" i N * 

X N { 2 1 




DC 60 

I r 2 , N 




IF (X.N( I )-XMjN) 

5?, 54 

, 3 A 

52 

XM1N e 

XN { I ! 




GO TO 

55 



54 

IF (X':( I J-XMAX ) 

5 C ,36 

,56 

56 

X'MAX 

X N "t I ) 



SC 

CONT I !\ 

ML' 



60 

COOT IN 

, [ *- 
’••• V 




P.ANG F. 

= XV AX- XV, 

IN • 



XCL = 

VC L 


. / 


WIDTH 

= RANGE /.< 

CL 



ACL - 

Cl.ASS+J .0 




XII) = 

>:.* '■ I N 




DC 70 

I = ? , '! C L 




I V 1 r 

3-1 




X ( I ) = 

X(l?i; 1 +v; 

IDT.-: 



X v ! 0 ( I 

't \ ) = x ( ! 

M 1 5 4 

I ATf 


7C COOT I ROE 

DC s: I = 2»-\Cl 
t y i *- i-i 

O ro t i( !-M J t C. 0 
I'O ?( J-- 3,N 

I f ( X *i ( J ) - X ( I 1 ) J 70 , 7? , 72 
7? irtv\( j }..;<{ ! , } 74,7:7,7s 


?;• 

IF 

( r \ ( 

» 

1 

- X( 

j 

) 

) 


7 4 

7 A 

Cf 

r- r '• ( 

1 

l 

• • 

1 ) 

t 


•s - 

1.1 

PL 

76 

cr 

NT IN 

■ J 

r 







00 

cc 

>:T 1 *; 


: 








or 

scot 

* ♦ 

V. 

L ) 

- 


t 

r. 

u 


DO 

u 3 

I 

r 

1 ,"C 

I 

i. 





i> 

?‘"r 

( 

l 

xm i 

r. 

( 

r 

• 

) 

-X 


if 

( : y 

p 

*• 

ico 

* 

) 


r» 

C» 

2 , 


Figure A-3 - 7072 Computer Program Listing for Failure Data 
Distribution Determination. 




Cl Lf KICm = v. 0 

PRINT 4-6. EXP.XMlCt 1) 

GO TO 33 
6? Cr-,v7 j 

CFRECU ) - AN*G. 3989 « 23 /X 51 G < ‘EXPi:FC“EXP)«v:]DTH 

63 CCNTUiCc 

PRINT 47 ♦ (XN( I ) . I - I »N) 

PRINT A?. » U FkLU t 1 ) »OFfcEG< I ) .X( 11 . h] .f'CL) 

M =•• MC L / ? 

J = 1 ' 

DO 84 I«]«K 

IF ( FFRF.C- ( ! ) - 4 . 0 » 77.64.84 
7 7 IF ( I -V ) 79. c 4. 64 
79 

I PI = 1 + 1 

t'FPEOdPJI = EFREOt I )+r.FR£(J( I p]l 
Of R£G< 1P1 ) = OF NFC ( I H-Cf'R5Q( jp} j 

64 CONI I NOE 
K = MCI 

mc tr = y . CL -' ■■;+ 1 

DC 87 hi .X 
MV I r. MCI.- 1 + 3 

IF CFF RFC (•'■•;: 1-4. C! 89.67.6? 

£9 JP (VMI-VfTf<) 37,67.86 
et X.e MCL-] 

MM 1 - W 1 - 1 

t" F RL C ( •*'.■•*. I ) * EFRECCMMI ) +EF.\fO ( 1 + ] 1 

CF RL O ( VV J i - OF RED ( i-*.M I > ,+OFR F 0 < V,M H D 
87 COM I MCE 

PP ! NT < t ?.t CLrKEOC I ) .Or RL'U( i ) .X( ! ) . I -J.X ) 


PR I 

NT 

4 8 

. ( 

MCI 

r 

K- 

j 

CHI 

SO 

= 

r> r, 

VI v 

DC* 

f.o 

i* 

J . K 

CM I 

so 

r- 

CHI 


CHINO t CHIF-c. ( (CJFKF'sJ ( I )-L'rp.CG< I 1 )*»CFRES( J J-LTFF C-» I 

j n/trPLOC) 

86 CCN 1 J.V.JC 

prvjRf: 

CAP 

? r oc r a v c c : t* f 1 7 ; i no Till m l a , si an d a r o deviation 

AND CO.'.r i i'.Ct l I PITS Or &AMXV VAPiAHl.CS f ROM 
ASSU-'ED .‘ivAVAL Di E Tr.' I CUT I CNN 
CO.N'F t 9 C. 

3 FORMAT (IS) 

9 r OP’-' A 7 inc.oi 
C FORMAT ((HFl 9 .cn . 

K fOP.-A: (Fie. ; , 13 / mf :o.on 

1? »C*«A1 flr'A// 10*5 SFECI-FNF, 7Z . *»•$ T *FJiG Tii . 1?X, 

; E'tU.OG '•*: AN. !«x . : lHL.OC. SI Df V , 7X ♦ £‘|F K.L'i.NCSN . 7>: . 


Pi «»ure A - 3 (ContM) 
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2 6HKURTCSiS/iXi!]0'F , 15.?»r2C'.&*)0H WJT.HN = * 

3 F?C.fi,2Fl5.5) 

n format ( i x » 2c.-iCo\ri3i,;ict Li.'- ns or » I 3 * 

1 J7H PERCENT ARE USED) 

1 A FORMAT (IX, i ?HCO.\F 1 DENCF. L 1 V l T S * 1 0 Fif-it AN I.CV/CR , 5X , 

1 lOHMEAN L'P^ERiSX* 13KSTO DEV LOWER » 6X , 

2 13MSTD DLV UPPE R/ 6X * F23 . » * F 1 6 . * , F 1 S . F , F ] 8 . *t ) 

15 FORMAT (A IX, ICHWJTH Nil ».F2 0.6) 

36 FORMAT (IX, 32K CHI SOLAR E FIT TC NORMAL DIST = , F14.6, 
1 5X , 2 3HD: GR.EE OF FREEDOM r ,1.3) * 

D 1 MEN £ J ON Xu (li'lMiX.M OG U 20 ) , CM It (100) , CM ! U ( 1 CO ) , T ( 1 00 ) 
PLAN 3.CCRF 


READ 6, (T( I) ,1=1,30) 

REAL 5, ARcA 
READ 6, (CM I LC I) , 1 = 1 , ICO) 
READ 6, (CHJ v< I ) * I = 1 » 100 ) 
16 RE/D '10, S,.N <XN< 1 ) » J--1 ,N) 
DO 2 0 i = i,n 
XMCGI I ) = LOGr (XN( I ) ) 

20 CONTINUE 
XNSl'M t C.O 
DO 2 5 1=1 , N 

X NAUM = XMSL.M 4 XNlGG(l) 

25 CONTINUE 
TN = M 
AN = M 

XMLOC = XNSUV/A.N 
XV2 = 0.0 
XM 3 =0.0 
XV.* = C.O 
DC 30 1 = 1, K 


XSO = (X.\L06( I l-XMlXC*)* (XNLCG* I l-XVLCG). 
XM 2 = XM 2-1 XSO 

XV?. * XH3 < XSO * ( XNL CG ( 1 1 -XMLOC, ) 

XV. A r- XMA-I XSC'XSC 
30 CONTINUE 

XS I G = SORT V ( XV.2 / (AN- 1.0)) 

XV. 2 r XV 2 /AN 
XV 3 = XV 3 /HI. 

XV A = XV 4 /AN 
XSiGL = SORIFU.V?) 

XX? = XV? /(XSIGL* XV? J 
XK4 - X V A / ( X.42 * .'<*•• 2 ) 

A = SC^T F ( AN ) 

IF ( .‘.-30 ) 3 6 , i 'i . 3 6 
3A x.Vl = XMLOC- - 1 (\-l ) -XSIGL// 

XVU - XVLOC-, 1 1 )» XSIGL /A 

GO 1C- 3f 

3 E X Ml = XV L OC- - A L A *• x S I G l. / A 


Figure A- 3 (Cont'd) 
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Xi 'U = XVLC6+ ARE A* x S ! GL / A 



Xi‘U = X 


IF (\-l 

35 

= IOC 


f-K - U 

37 

CCU 7 ! Mi 

35 

XT-TPL = 


35 XT-TPL = X5lCL*SC-MTr((AN-].C)/CH!L(N-l) J 
X5 7PU = XSiuL*SGRlff (AN-I.OJ/CMIUCN-I ) J 
N = lfi 

prim i?i n»s*x;-:log»xs!Gl* x>;3»xx4 

PR1KI 15» XilG 
PRINT 13 » COf-.F 

PP INT J4, X.-’L , X MU , Xc‘ I PL , XS I DU 
IF ( N- 1 0 ) 40.4C.39 



39 CONTINUE 


call >:in 


PRINT 16. 


40 CONTINUE 

\ 

GO TO IP, 

i 

END 


i » '.w u 

CAL L < I F I T ( ft » XKLCG* X K 1.06* XS I GL .CHI SO.M C 


I 


Figure A-3 fCont'd) 


APPENDIX B 


FORTRAN Computer Program To Determine Time Dependent 
Strength Distribution parameters 


B-l . Flow Chart 

B-2. Definition of Program Variables 
B-3. Computer Program Listing 




Calculation of failure 
distribution parameters 
for interpolated stress 
levels 


Figure B-l 


Flow Chart for Time Dependent Strength Distribution Generation 
from Failure Distribution Parameters. 












I 





r 


Cummulative lognormal 
failure distribution at 
N=LK stress levels 
from failure distributio 
parameters . 


Figure B.l (Cont'd). 











Strength frequency 
distributions for 
various log cycles 
of life. 


Parameter calculations 
for normal distribution 
fit to strength 
distributions . 


Figure B-l (Cont'd) 
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M = TOTAL NUMBER OF GENERATED STRENGTH 01 ST. 

' IfiT = NO. OF VALUES OMITTED BETWEEN PRINTOUT 

XJMIN -• MINIMUM. LOG CYCLE VALUE FOR STRENGTH DlST. 
XJINT e INTERVALS BETWEEN LOG CYCLE VALUES FOR 
STRENGTH DISTRIBUTION S 
N = NUMBER CF FAILURE DISTRIBUTIONS 

SINT r. MINIMUM STRESS INTERVAL BETWEEN INTERPOLATED 
FAILURE DISTRIBUTIONS 

A STS ( J J = STRESS LEVEL FOR EXPERIMENTAL FAI LURE 
DISTRIBUTION 

AM LOG ( J ) ~ EXPERIMENTAL FAILURE* DISTRIBUTION MEAN IN 
LOG CYCLES AT STRESS AST R 
ASIGL(J) * EXPERIMENTAL FAILURE DISTRIBUTION 
‘ STANDARD DEVIATION IN LOG CYCLES 
DSTS = DIFFERENCE BETWEEN Tv/O CONSECUTIVE 
EXPERIMENTAL STRESS LEVELS 

DMEAN - DIFFERENCE BETWEEN TWO CONSECTIVE EXPERIMENTAL 
FAILURE DISTRIBUTION MEANS 

. OS 1 G = DIFFERENCE BETWEEN TWO CONSECTIVE EXPERIMENTAL 
FAILURE DIST. STANDARD DEVIATIONS 
SLAM(J) = SLOPE BETWEEN TWO CONSECUT I VE EXPERIMENTAL 
- . FAILURE DISTRIBUTION MEANS 
SLSIG(J) = SLOPE BETWEEN TWO CONSECUTIVE EXPERIMENTAL 
FAILURE DISTRIBUTION STANDARD DEVIATIONS 
L = COUNTER FOR INTERPOLATED FAILURE DISTRIBUTIONS 
MX = DSTR/SINT 

LX = sequential numbering of sample failure 

PARAMETERS AND INTERPOLATED PARAMETERS 
STR = STRESS LEVEL FOR I TH FAILURE DIST. PARAMETERS 
XMLOG(I) = FAILURE MEAN IN LOG C,YCLES FOR 
EXPERIMENTAL AND INTERPOLATED VALUES 


XSIGL(I) = FAILURE STANDARD DEVIATION IN LOG CYCLES 
FOR SAMPLE AND INTERPOLATED VALUES 
K = L X = NUMBER OF GENERATED FAILURE D I STR I BUT 3 ONS 
BY interpolation 

XVAL e LOG CYCLE VALUES FROM XJMIN TO XJMIN ■» M*X J 1 NT 


FOR M STRENGTH DISTRIBUTIONS 
2 = STANDARDIZED DISTANCE FOR LOG CYCLE FAILURE 

DISTRIBUTION AT I TH STRESS LEVEL FOR K TH LOG CYCLES 
PRC; 3 = NORMALIZED PROBABILITY DENSITY UNDER LOGNORMAL 
D I S T R 1 B UT I ON FROM - 2 TO 2 

AREA (NX) = CUMM’JLAT I VE NORMALIZED AREA OF LOGNORMAL 


FAILURE DISTRIBUTION FROM ZERO TO 2 
TSTRMNOJ "STRENGTH LEVELS PRINTED OUT 
HINT UNO) - IM EAGER REPRESENTATION OF TSTRUNO) 
CY C(<) r LOo-CyCLE VALUE CORRESPONDING TO K S 


Figure B-2. Definition of Program Vatiablcs. 
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FREQ { 1 i K ) = FAILURE PROBABILITY FROM CUMMULAT1VE LOG 
NORMAL FAILURE D I STR I BUT I ON AT LOG CYCLES CYC 
BE 1 WEEN STRENGTHS 1(A) AND 1(B) 

UJ = ABSCISSA VALUE OF STRENGTH HISTOGRAM CLASS 
INTERVAL 

FI - FREQUENCY OF OCCURRENCE IN STRENGTH HISTOGRAM 
CLASS INTERVAL 

F I U3 = FI TIMES UI 

SMLCG ( .< ) = MEAN STRENGTH AT SET LOG CYCLES 

SM2 i SM3 i SMA = 2ND. 3RD. AND ATM MOMENTS OF 
STRENGTH DISTRIBUTION 

SSlGLLK) = STANDARD DEVIATION FOR STRENGTH AT SET 
LOG CYCLES 

SK3(k) : = COEFFICIENT OF SKEWNESS FOR STRENGTH AT SET 
LOG CYCLES 

SKA ( K ) = COEFFICIENT OF KURTOSIS FOR STRENGTH AT SET 
LOG CYCLES 

TOTAL NUMBER OF I NO S - N/NINT WHERE N = LX MAX* 

TOTAL NUMBER OF J S = N SAMPLE FAILURE DIST. 

TOTAL NUMBER OF I S = N WHERE N = LX MAX. 

TOTAL NUMBER OF K $ = M 


Figure B-2. (Cont'd). 


o a r> 


183 


*** SN! 1 T H * f E 

4 c p Ilf r c. ?. r s a \ , e x e c u t ti r o r i r a n " 

SUDS' CUT I i\l PRO- (2 .PRO L' ) 

c P.-C'3 * THl area l-holr normal DU.tr I CUT ion BETWEEN 
C PLOS AND MINUS 2 STANDARD DEV i AT IONS 

IF <2-1. 2) 1 000 » 1 000 » 10 1 0 

lore 2 sc * 2*2 

. PK03 - 0.79766^6*2* ( 0. 9599977A-2SG* (C . T A659A 33 
1 *-2Sv7*{ 0.02463C3 10 - 2$0*C , 90? 357 ACS7 ) ) J ~ 

GO 7 C 107.0 

1 C> . I F ( 2 — ? • a ) 1 w 2 1' * 1 0 0 0 » J 0 6 0 
1 020 2 SO - 2 s 2 
PFOB - 1.0 
PTERM =1.0 
fact =1.0 
ODD I h ~ 3,0 

1020 PTE ■•' = -PI Ek.‘.*2S0/ < 2. 0"F AC T 1 

T * F.’.'l = PT r. Rf 'i/OODJ \ 

PROS = PROD -4 TERN; 

, * F CA^SFJTtPMl - 0.00007) 1 090 .1040*1 GAO 
10-4.- r ACT = r ACT + i.o 

Ot.-;.* I A = COD IN + 2.0 
CC: TO 

1000 PROD ' = £.797asV*5*2*PRcp' 

GO TO 3070 
1060 PEC = !,/(?_< 2) 

FPClI = 1 . -0. 75738AS'i*tXPEF (-2*2/2 .0) /2« 

' a-* V R£C * U * - r EC*( 3. - RtC* (15. -RE C *105 . ) ) ) ) 

lu.'v t I l.'.t'c. 

Return • • 

end 


K - 23 

’ AMLOGdOl. AS I G L ( 10, » SLAM (9) * 

i oi- ( 253 > » APEAC2S3 ,235 , FRED ( 2 K . 1 . ?V 

3 ( ’ ^:CL(23»t S.<3f23), SK4(?3)* 

j? *' r. ; : * 1 V - x 1 1 -■ J 1 J *' C. ( 1 0 , 2 3 ) , GFR £ 2 ( 1 A , ? 3 , 

CF — •• 

? <wU*CFrtSQ)!\I^ 

5 READZO , x.w I N i >; J I t ’ * U 

1C FORHATt ?F 1C .0) 

FAf !‘Au.rc? 5 ?Z ,? i BUTION VALUCs ARE READ lH FROM 

lowest iO highest stress and Stresses are 

IRUSER VALUES OF SINT 
Re AO 20 

riM,re 8 ' 5 - Du\^„ er 4 r :So„ Ustins for *"•** 


t~\ O f\ 


20 (50!-. JSLANK 

,, r c ’ j) ,an:loc-< j> .a$jsl(J) »j«i 

30 F Cn ••‘.A " i ( ! & » r £. • C / { 3F 1 0 • 0 )) 
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DO !8C J=1,V. 
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CO 1C 130 
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1 I3.F9. 6* 

I 2 » 

F7.6, 

i 3 » F 9 . £ , 

I 3 . f 

r 9. 

6*13. 

F 9 , 6 ) 


STP = STR+ 

5 1 \ 

1 * X l h 

T . 





2 2 0 

CC.M 1NUE 









224 
22 2 


DO 22? J-l.NOLD 
1 * INDEX (J) 

x.\i?» = oor-itJ) 

DO 22 A K*l,y. 

OF RE (M J , £ » - x ■ ,U-1 < A ?< A ( I , \ ) 

comi-ive 

ccm - • 

C. fPEOCFNCY DiSl A.':?. NORMAL [MSTRirTJT ICN para 
XJVl.'i = XJ/IN-XJIM 
CO 260 Jt’.:-; 
xj = J 

CYCU) - X J fc X J ! !\ T * X J*’ I ;; 

PPiN] C . J. CYC(J) 

2 AC 1 f r/<V.Al ( JkT.-V-.J = *13 » 5X » 1 3 ML Go CYCLtS - . 

PRIM 250* (1 tAVLAd • J) •!>] IM1 
250 FCSXA r ( 2/ . jti ( I i . 1 /. , 9M? :<l vil'f \C Y / i 1 X , l . , f 1 0. 


I C- . f , , 1 5 


G, !'.*•? K .6, !5.f 1 


c- » 


1 1 


260 CON! I AL E 


0 V.EAf; » 

= * F 5 ♦ 6 ) 

AT 1 VL , 

9 . 6 , 


v t'TERS 
Ft. 3 ) 

5. ! 5, F 10.6 
1 0 . 6 . ! 5 , F 1 0 . 6 


Figure B-3 (Cont'd) 


NVi - ::-i 
DO 270 JO ,!•' 

CO 2 7 0 ! » NV.l 

ir*i = n i 

FF.EOSiiJJ = A T. il A ( I ? 1 » J ) -AREA < I » J I 

?7cj ccrtinuf 

m = ,\”1 

NO.nM-L 0 I STK I CUT I or; P ARAKS. T FRS OF HIST OGRAV 


ML A a A AO 
DO 3 CO J*] 


FJ * o . 

r i u i - o . 

UI = SM!M-sn;l*C .5 
DO ?B 0 1 = 1 ,M 
UI *• Ul+SIMT 
FI = FUFRF j( I * J) 

M UI * f I U I +F RLO ( I * J ) *U I 
2 SC COAT I RUL 

SMLOC(J) - FI'JI/FI 
SM2 = 0, 

SV3 =0. 

SMA =0.0 


UI = 

S 1 .\ - S I N T * 0 • 5 

DO ? 

90 1 = 1 »K 

•J5 = 

uh-sim 

SC = 

tu: -facet Jj i 

V2 

= S*-'2 4 SO* FRF 


* S" 3 4 r'KfU ( I * J ) « ( U J -Sf'LCG ( J ) I * SO 
SVA r S v .A4S«‘*S0*FffLQ{ I ♦ JJ 

?vC COAT I A'.;* ‘ 


S."-2 

' r 

SM2/FI 








SAT 3 

, r 

SV.3/F 

I 








S.v.A 

r 

S'-'A/r 

I 

• 







SSIGKJ) = 

SORT f ( G."2 ) 








SX 3 

( J5 

r SM 

13/ISSIGl (JJ* 

' Sr 

:2 » 






' SKA 

( J) 

= S> v 

'4/(S:-:2*SM2) 








cc\ 

T : 

! : ' r 








, • 

-P I 

•\ T 

?? 







• 


Pn I 

AT 

210 






. 



•? ff.iV 

■ .•% i 

{ 1 Of:! 

.\U-;-.‘:K .1 

Oil co 

c y 

CL 

t'S . sx 

i 


1 

i fr: 

T. A' ; ST K l o T H » “• 

x» 

I SH 

r t 
t 

^ « 

DCVI 

Alio: 

N • ! OX . 

2 

3 Or! 

5.KF 

.v\LSS . j OX , 10 

H 

•<!J 

P T 0 

J 4 

0 ) 



PF 1 

* ‘ T 

• % . 

•20. 

u. crct j) , 

S* v 

LOO 

( J 1 

« 

SS I C-L ( J 1 , 

S.<3 ( J J 

j 


{ Jl * 

o r ; . ) 







PF:J 
c* rvR 

.Ml 

y. a t 

32 0. 
((Ir 

( J. CrC ( J) , 

» r 9 « 2 . 2 ? ? A * / 

.2 

l OS 
F?0 

( jy 

. F:> 

t 

r 

SSIGl. 

( J ) . 

SK3 C J) 

pr i 

M 

32S . t 

AU-.J j; , I APEX 

( J 

) . J 

- j » 

‘iOL D ) 



s TCP 

« 

V a T 

p* /• I 

( •' 

x. 1 S.-i.MiSfRv 

j D 

M* 

-/ .F F 

A . 

?X» 

n hi; 

IT EOF R 

I 

1 1 : 

X » 2 1 2 

•o ) ) 







f» m 

see 

K*- I • 










Figure B-3 (Cont’d) 


4 { J) 




! i 


! T 
! 


it i 

i 4 ! 


187 


ic 4 00 J= j » .NOLO 
X * INDEX (J) 

S T ' J ) = SM ! .*\+ S i NT * X -$ I U 7 

2 = (ST ( J}-£ v.LOG (<) J/SSIGL K) 
1 f (Z) 320.3^0,3^0 
330 2 V -2 

I ? c 2 - . r> } 3A 0 , 3*i 0 1 360 

340 CALL P^03I2.PRD5-) 

PRC BA - C . 5-PR03-* C . 5 
GO TO ?. ?C 


35C 


I f 7 ( 2-3 . 3 } 35.5.2 53,365 


CALL PTDru.PP.Cc.) 

P ROB A = . P R.C6 « 0 . 5 + 0 . 5 
GO TO 370 
360 PRC BA =0.0 
GO" TO 3?C 
^c^.* r\-.C3A — 1 « C 

3 70 GF>£'6{J,.<| = PRCb A 

X':UM = MIJM(J) 

EFPEGIJ.JC) = G.-.7E:.-{ J,K)*x?:uv 
CMISOC Jf r. ) . = (OrRE-j{ J,K)~EF*?EC{ J,.< ) )*(OFP r C( J 
1 -FFRtC ( J,K ) ) /£ FKEQI J,K) 

*ec- C OMT I t\;:jE 

P R I I /t 0 0 , K 

4 90 FOR' AT ( 1 HO , 3H K=,J 5 ) 

, ?R i S 1 ^’ tST( ' ,, ’ CF ^ C,J,K, » 0FRW(J ‘^’^ RFC { . 
J ».f<ISQCo,K} , J=1 ,f;0L0) 

A^5 FC.?i>;AT ( 16X, i AhSTAiiftGTIi LEVEL, 9X , 1 ] tiNOR"-*Al FRJ 
; 7X,13HLXPFCTL0 FREQ, AX, 

SC-o'c-'vn.v-Vr ^ UAR " VALUL/ * 10X»KC20.'6 ) ) 

GC TO 5 " 


,KJ 


1 . y 


K ) , 


7X 


r K 


-.i-- 


Figurv B-3, . (Cont’dJ. 


188 


APPENDIX C 

Listing of Fortran Computer Program CYTOFR Used To Determine 
Cycles To Failure Distributions 


C-l. Flow Chart 

C-2. Definition of Variables 

C-3. Fortran Program Listing 


Program to Calculate Parameter Estimates for the 
Normal and Log-Normal Distributions and Conduct 
Goodness-of-Fit Tests 

MAIN PROGRAM (CYTOFR ) 



Fig. C-l 


1 

















m 



i 

i 












Subroutine to Find Area Under Standard Normal Curve 

FUNCTION PR03(X) 


IS (X-1.2) 


XSQ = (X) 2 


PROS = ( .79788455 ) (x) . 99999??4-XSQ[ .16659433- 
XSQ( . 02463831 0~XSQ( . 002397486?) } ] 


0 . 


is x-2.9 


RETURN 


XSQ ~ X* 

PROS =1.0 
PTERM =1.0 
FACTOR =1.0 
ODDI I IT =3.0 


RECXSQ = -~ 

PROS = 1 .0-( .79788453) EXP[ 

y2/ p 

) ( i , 0 -RECXSQ ( 1 . -RECXSQ 
( 3 . -RECXSQ ( 1 5 . -RECXSQ (105)))))] 


0 

^RETURN 

C -1 (continued) 
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Fig. C—1 (continued) 
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Subroutine to Find the Moment Coefficients of 
Skevrness An<3 Kurtosis 

SUBROUTINE ALPHA 


Fig. 






List of Definitions for Program to Pit Normal 
and Log-Normal Distributions to Cycles- 
to-Failure Data (PROGRAM CYTOFR ) 

Main Program j 


NDATA DATA 

= number of observations. 

STRLV 

stress level in psi. 

AKURCY 

accuracy to vrhich cycle s-to-f allure data 
are knoun. 

RATIO a 

stress ratio 

X(I) 

cycles-to-failure data 

CUMFRQ ( I ) = 

cumulative frequency of each X ( I ) j ie, 
number of X's less than or equal to X(I). 

PCAPtEA ( I ) -- 

C U HFRQ ( I ) /NDATA 


Subroutine to calculate the mean and standard deviation of 
the cycles-to-failure data (SUBROUTINE MEAN) 

SIGMA = sum -of the X(l)'s 

XMEAN s= average of the X ( I ) ' s 

n 9 

T0P2 « £ (X(I)-XKEAN) 

i=i 

DEV = standard deviation of the X(l)'s 

Function subroutine to find the area under the normal 
curve (FUNCTION PROB(X)). 

X s abscissa value for vrhich corresponding area 

is desired, 

PROU = desired area,. 


C-2. Definition of Variables 
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Subroutine 

for 

Chi-square goodness-of ~f it test (SUB- 

ROUTINE CKISQA). 

K 

8 

number cells. 

XMAX 

8 

largest value of cycles-to-fallure . 

XMIN 

8 

smallest value of cycles-to-failure , 

CSV 

ss 

cell starting value. v , 

CEV 

ss 

cell end value. 

CLB 

S3 

cell lower bound. 

CUB 

8 

cell upper bound. 

FREQ(J) 

8 

J.1. 

number of observations in J cell. 

BEQAREA(J)^ 

j-y. 

expected value of J cell. 

CHISQR 

-- 

total Chi-square value. 

u(i) 

8 

th 

Chi-square value of I cell. 

Subroutine 

for 

Kolmogorov-Smirnov test (SUBROUTINE DTEST). 

z(r) 

8 

abscissa value on standard normal curve for 
a given X(I ) , 

ABUNCN 


area under standard normal curve from - 
to Z(I). 

DSTAT(I) *= 

absolute difference betvreen the data cumu- 
lative frequency and the hypothesized cumu- 
lative frequency. 

. XHEAN 

8 

average of the X(I) 's. 

DEV 

S 

standard deviation of the X(I)*s 

PROB(T) 

= 

area under the standard normal curve from 
-T to +T. 


C-2. (Cont'd). 


Subroutine to calculate the moment coefficients of skew- 
ness and kurtosls (SUBROUTINE ALPHA), 


ALPHAS 

- 

moment coefficient of skevmess. 

ALPHAS 

SI 

moment coefficient of skevmess, 

VAR 


T. (X(I)-X) 2 



1 

xH 

If 

•H 

T0P3 


E (X(I)-X) 3 



1=1 

SKEW 

= 

third moment of the data. 

STDEV 

n 

biased estimator for standard deviation 

TOP4 


n h 

E (X(I)~XF 



1=1 

TXURT 

— 

fourth moment -'of the data. 


C-2. (Coiit 'd) . 


I. 

f 0003 


.• ..003 

! o' o o 3 


\ 0 0 0 7 

I 

00025 

; 00?5 

0030 

00043 


00043 

00^5 

00 4 7 
COO 50 
S0051 
0 0 53 
000^4 
000 57 

go 6 c 


0 0 0 63 




DC 65 
070 

v» 0 7 1 


CO 07 3 



00101 

*010? 

0106 

oon? 


= 011 ? 
,>0113 
00l?3 

>01 ?3 
P0)?4 
0 013 ? 

0013? 
CM3 6 
'^136 
'0 151 

151 


PROGRAM CYTOFri ( INPUT i OUTPUT » T APE 1 = 1 NPUT) 

C- PROGRAM TO FIT NORMAL AND LOG-NORMAL CUR\'E TO DATA AND CHLCK 

C GOODNESS OF FIT. 

DIMENSION X<100> » CSV ( 9 ) , CEV (9) iCLB(9) iCUU(9> » CUMFRQ ( 1 00 1 i 
IPCaREAUCO) »DSTaT(100) »KRE0(9> «APEA (9) » REDaREA l 9 ) » EXFRlQ ( y) »U(9> * 
22 ( ICO) . N X ( 1 0 0 ) iRAMK(IOO) 

EXTERNAL prod 
7 1 o PRINT 1 

c — - mdata=data=number of observations 

C STRLV = STRESS LEVEL IN PSI. 

READ 6 »MDATA,DATAiAKURCY,STRLEV, RATIO 

c x= number of cycles to failure 

6 FORMAT (I 3 »F 5 . 1 »F 9 . 4 »F 10 . 1 » F 8 * 5 ) 

IF (EOF i i ) 56 i ^5 

55 READ 7 » ( X(I) i I- 1 iNDATA) 

7 FORMAT ( 6 F 1 0 - 0 ) 

C 

C SORT X ( I > TERMS in ASCENDING ORDER. . 

C 

K-NDATA -1 

IF(K.LE.O) GO TO 30 
DO 20 1 = 1 ,K • 
n=nqat a~i 
1 STOP =0 
DO lo J=i»M 

IF ( X ( J ) • LE • X ( J+ 1 ) ) Go TO lo 
SAVE=X(J) 

X(J)=X(J+ 1 ) 

X ( J* 1 ) =S/i.VE 
ISTOP=ISTOP+l 
1 C CONTINUE 

IF (I STOP. ED. o) GO TO 30 
20 CONTINUE 

C. 

C SET CUMFRO(l) ARRAY 

C 

30 DO 40 I=i,NDATA 
40 CUMFRQ(I) = I 

C PCAREA e F (N) OF OBSERVATIONS 

00 759 I = ), MOAT A 
75 g PCAREA ( 1 ) = CUMF R 9 (I ) /DATA 

Pr<r j T 405 

4 C 5 FORMAT ( 4 c*» 53 HN 0 RMAL DISTRIBUTION FITTED TO CYCLES-TO-FA 1 LURE D^ 

lA.///) 

IF (RATIO. eQ.C.O) GO yO 414 
PRINT 4 0 ?, STRLEV, RATIO 

402 • FORMAT (? 9 X, 1 AH STRESS LE VEL= , F 1 o . 1 i 5 H PSI. » 16 X 
i,] 4 hstress Ratio -,f 6 . 3 //j 

GO yO 415 

414 PRINT 416, STRLEV 

416 FORMAT (? 9 X, I 4 HSTRESS LEVEL* * F 1 0 . 1 « 5 H PSI.»l 6 X 

i,? 3 hsiress Ratio = infinity//) 

415 PRI’lT 4 C 4 

4 C 4 FORMAT ( 55 X, 22 HCYCLES TO FAILURE DATA/) 

PRINT 40 ?, <X(I) , 1 = 1 iNDATA) 

4 03 FORMAT ( 6 (] 0 X » F 1 5 , 3 } ) 

PRINT 3 


C-3. Fortran Progam Listing. 






3 


00Q155 

{ : 0)56 

I "' 016) 

000171 
r ^0201 
>o ? i c 

OOQ212 
0 2 1 3 

.\,?2 6 
&V023C 

0 C 6 2 3 4 
ITio?4 0 
f fi0?4 0 
01)5244 

1 ) 0 ? 4 4 
100245 
CC 0255 
| )Q256 
I J0 ? 64 
000270 

' I I1027C. 

l!’0303 

000303 

S > 0 3 0 7 
:o3i3 
/ 0323 
C C 5 3 3 3 

t >p343 
;o343 
000345 

i 


I 



54 


1 

401 


417 

418 
2 

4 1 3 


56 


FORMAT ( 1 H 0 ) 

CALL MEA”(Xi Data* NDaTAi XMEaN» DEV) 

CALL CHISQa<X* 1>aTA» NDATa. PROb* AKURCY* XMEAN. dev, 2) 
CALL D T F. 5 T (PCARE 4 , NOATA , x « DEV * DST AT , PR°U i X^EAN , 7 .) 

CALL ALPHA (X, NDATA, DATA, XMEAN, DEV, ALPHA3, ALPMA4 ) 
AKURCY = * 0 C 0 0 1 

DO 54 1=1, NDATA 

NX(I) a(ALOG(X(l)/2O.)*AlOG<20. >>*100000. + .5 

X(I) « NXiI> 

x(i) = xti>/ioo66o. 

PRINT 1 
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FORMAT <lHi ,3HASD/// ) 
PRINT 401 

FORMAT ( 30X » 57 MI.OG- NORMAL 


DISTRIBUTION FITTEO TO CYCLES- TO-FAILURE 


} DATA.///) 

IF (RATIO. EO. 0.6) GO TO 417 
PRINT 402, STRLEVj RATI 0 
GO TO 410 


PRINT 416, STRLF.V' 

PRINT 2 

FORMAT (49X » 34HL06S OF THE CYCLES TO FAILURE DATA/) 

PRINT 413, ( X C I ) ,1 = 1, NDATA) 

FORMAT ( 6 { 8X , F 12 « 5) ) 

PRINT 3 

CALL MEAN (X i DATA, NDATA, XMEAN, DEV) 

CALL CHISQA<X, DATA* NDATa, PROB» AKURCY, XMEAN, DEV, 2) 
CALL DTEST (PCAREA, NDATA, X»DEV»DSTAT,PR°B, XMEAN, 2) 

CALL ALPHA ( X , NDATA, DATA, XMEAN, DEV, ALPMA3 , ALPHA4) 

GO TO 710 . 

STOP 


END 


nimr m*m* • ' «* • 


SUBROUTINE MEAN (X, DATA, NDATA, XMEAN* DEV) 

C --SUBROUTINE to calculate the mean and standard deviation of data. 


f ’>0010 


DIMENSION X(NDATA) 


« 0010 


SIGMA* o.o 


000011 


DO 8 1=1* NDATA 


i \6oi? 

8 

sigmas$igma* x ( i ) 


{ 016 


xNEaM s sigma/data 


Co-0017 


TOPS = C.O 


OOO0S0 


DO 9 1=1* NDATA 


1 ,’0021 

9 

T0P2 = TOP2 ♦ { X «: I > - 

XM£aN) 

C 0 C 0 2 6 


DEV =5 q RT (T0P2/ (DATA 

- 1.0)1 

000036 


PRINT 14, XMEAN 


d ]: 5 0 ^ 3 


PRINT 15. oEV 


d ;>oo5a 

14 

FORMAT ( 1 OX * ISHSANPlE MEAN=» F17.6) 

000054 

15 

FORMAT ( 1QX» 15HSTD. 

DEVIATIONS* F14.6) 

<1 ‘0054 


return 


fl >0055 


emd 



Is 

l 

I 

I 

E 

I 

II 




C-3. (Cont'd). 


— — - aggagisaaBBiMr 




nooo3 

4 |Q005 
t . 006 


i i> 0 0 1 6 

6^0017 
000022 
fi 1*0023 
(| 1*0024 
000025 
0*0026 
4 *0030 

Ov-0034 

000035 
< )0037 
i 0042 
OC$044 
4 10046 
<J s 0046 
000051 
Qpj>051 
d I* 0 0 5 3 


000075 
<f O0? '6 


FUNCTION PROU(X) 

C THIS SUBROUTINE GIVES AREA UNDER NORMAL CURVE FROM -2 TO +2 

C WITH AN ACCURACY OF 0.00005 208 

C ~-Z VALUE GIVEN BY CALLING PROGRAM MUST BE A POSITIVE NUMBER. 

IF (X-2»?) lltllf 12 

11 XSQ=X*X 

PROB= 0. 797884550X0 (0.99999774-XSO«Mc.l6659433-XSU<M0.024b3S310-X 
1Q«0. 0023974867) ) ) 

RETURN 

12 IF (X-2.9) 13,14*14 

13 XSQsX*X 

PR08*=1 , 0 
PTERMsl.C 
FACTORsl.O 
ODDIHT »3 . 0 

970 PTERM=“PTERM«XSQ/(2.0«FACTOR) 

TER m =PTER m /°OOINt 

prob=prob+term 

IF ( ABS (TERM) - 0,00007 > 60,90,90 
90 FACTOR *FacTOR*1,6 

0001NT=ODOINT42,6 
GO TO 97c 

8o PR0B = C*79788455«X«*PR0B 

return 

14 RECXSR= 1.0/ (X*X> 

PROOs 1.0 - 0»79788453«’EXp(~XoX/2.C)/X*(1.0"RECXSO#(l, -RECXSQ*<3 
1 - RECXS0<M15. - RECXSQ# 1 05. >>)> 

RETURN 

E n O 



f 


: 1 


C-3. (Cont'd). 


i 


L 


0 0 0 C 1 3 
: 013 


000013 

<ji00?4 

000025 
0^0027 
(j 0030 
Fo6o32 
000037 
(j 0046 


C j 0050 

c 6051 

000057 
6*0060 
d 0062 
Ou0067 
o 6 jo 75 
dlol 03 

cL .51 10 
conin' 
Ojiol 12 
46117 
0T0123 
C'l27 
ffv-135 
W0l37 

o66n3 

J 0144 

6^50 
000152 
Pi 53 
0172 

Oi6l77 
011620? 
0»i5206 
C00213 
Cf 0220 
01 0?24 
000243 
000244 
0 0246 
0 v0?47 
000250 

o' ??7? 
CCO300 

Of 5300 
‘I ?3 1 (1 


SUBROUTINE CHISOA (X, DATA 1 NOATA, PROD* AKURCY * XMEAN, DEV,/!) 

c SUBROUTINE To FIT a HISTOGRAM TO THE data ANO PERFORM the CHI-SO'JaK£ 

C-- TEST FOR THE NORMAL OR LOG-NORMAL DISTRIBUTIONS* 

DIMENSION X(NOATA) »Z ( f JD A T A ) ,CSV(9> »CEV(9) »CUM9> ,CUB(9) t 209 

IRE CARE A (9) * AREA (9) iEXFREQ (9) » FREQ(9)» U(9) 

CHlSORs ,0 

C— TO DETERMINE THE NUMBER OF CLASS INTERVALS, K 

K= 1.5 ♦3.322°ALOGlO (DATA) 

REALK=K 

c in order to determine the range, find x(max) and x(HIN) 

XMAXsX(l) 

XMlNa X(l) 

DO 17 1*3 , NDaT A 

IF C X(I).GT*XMaX ) XM A X s XU) 

17 IF (X ( 1 ) .LT , XMIM) XMJNsXU) 

RANGE= XMAX- XMIN 

C- TO DETERMINE the CLASS INTERVAL WIDTH.W 

c — --routine to round off class width to same number of places as the acc 

DIVIDE s 1 . 0/AKURC Y * 

KW = (( <RaNGE*AKURCY)/REALKN,5*AKURCY)*0IVI0E 

RK1 = KW 

W *= RK1/DIVIDE 

PRINT 62, XMAX 

PRINT 63, XMIN 

PRINT 65, W 

DO 22 1=1 ,K 

A= I 

B = C.5«AKURCY 
CSV C 1 ) = XMIN* (A-i.O)*W 
CEV(I)- CSV (I) 4W-AKURCY 
CLO ( I ) = CSvU)-B 

22 CULM I ) = CEV ( I ) *B 
CEV(K) e XMAX 
CUB(K) = CEV(K) *D 
0° 23 J=I,K 

23 FREQ ( J) =C • 0 

DO 24 I cl f ND aTA 
DO 24 Js j , K 

IF ( X(I) .GE.Cl.B(J) .AND. X (1 ) ,LE ,CUB ( J) ) FRECM J) *FREQ ( J) * 1.0 
24" CONTINUE 

c chi-square test 

PRINT 41 
PRINT 406 
DO 30 I = l'»K 

Z<I)=( CUB ( I ) - XMEAN) / OEV 
T = ABS( ?<I> ) 

30 AREA (I ) = PROB(T)/2.0 
REOaREA(I’) = 0.5 - AREA(l) . 

manu=k-) 

DO 32 1*2, M ANU 
M— 1 — 1 

IF( (z(l) ,GE, 0.6, and. 7 .(M) ,GE, 0,0). Or, ( * (I ) .LE , C . 0 , AND , 2 (M > # l-E , n . 

1 ) . ) GO TO 31 
REQAREAU)*: aREaU)*apea(m) 

GO TO 32 

31 REQAREA(l) = ABS < AREA ( I ) -AREA (M) ) 

32 CONTINUE 



C-3. (Cont'd). 



«r v <. « - a w« / v v V * • \ c- * ' % ♦ » 4 < 

0 0 0 317 DO 86 Ms 1 , K 

(T 0320 EXFHEU <M) sDAT A*RE0ARE A (Ml 

fl 6323 U(M)s(( EXFPEQ(M)-FREO(M) > «*2) /EXFREQ (M) 

000331 80 CHISOKsCH1SOR*U(M) 210 

c to print the table for chi-square test 

(j 0337 DO 33 I-lfK 

" T341 33 PRINT 34,1, CLR ( I ) tC^B < I ) * EXFREO ( I ) ♦ FREQ U ) * U ( I ) 

l>v04 1 2 PRINT 35, CHISOR 

q 0420 62 FORMAT ( lOXi 1 4HMAX I MUM VALUE= i F 15. 6) 

q 0420 63 FORMAT ( loX* 14HMINI MUM VaLUE=» F1S.6) 

000420 65 FORMAT ( 1 OX * 1 2HCLASS WlHTH=, F17.6) 

6° 0420 41 FORMAT UH0) 

C 0420 406 FORMAT <eX,5M CELL , 1 OX , 1 CHLO«E« CELL , 1 1 X , 1 OHuPPLR CELL » 1 3X » 8HEXPEC 


' 1TED , 13X» 8H0BSF.RVED» 1 3 X * 1 i HC MI -SQUARE D/8X , 6HNUMBER , 1 OX » 8H BOUNDARY » 

213X, 8H80UNDaRY,13X»9HFREQUENCYi12X»9HFREQUENCY»12X,13HVALUE OF CE 
. . 3LL/) 

r 0420 34 FORMAT ( 1 C-X , 12 1 5F? 1 , 6) 

006420 35 FORMAT Uhq*81Xi25hT 0TAL CHI-SQUARED VALUE =tF10.6) 

C 0420 RETURN 

d 0421 END 

I. 

[ 

I 

t 

[ 1 
I F 


I 

[ 

fi 

|| C-3. (Cont'd). 


' SUBROUTINE DTEST { PCAREA , NOAT A , X i DE V , DST AT , P- ' ,XMEAN,Z> 



r ■ 

c — 

— SUiiROUTIiiE T 0 CALCULATE THE KOLMOGOROV-SMIRNOV D-VaLUeS, ■ 

0 

3 013 


DIMENSION PCAREA(NDATA) »X(NDATA) ,Z (NOA1A) , DST AT (NDATA) 211 

c 

.10013 


00 706 1=1, NDATA 

f/( 

10014 


2U) = ( X ( I ) - XMEA n )/DEv 

0 

■) 0 2 1 


IF (Z ( I ) ) 703, 704, 705 

f 

0024 

7q3 

T =ABS (2(1)) 



C--- 

— ARNCNsAREA UNDER THE NORMAL CURVE TO LEFT OF Z FOR NEGATIVE 2. 

*\ 

111.027 


ARUNCN = (l.C-PROH(T) )/2.0 

0 \ 

J04 3 


DSTAT(I) = ARUNCN - PCAREA(I) 

00 0050 


GO T 0 706 

0! 

p>050 

704 

DSTAT(I) = .5 - PCAREA(I) 

0 

[1)054 


GO TO 706 

0< 

51)054 

7o5 

T = 2(1) 



c— - 

— ARUWCPsAPEA UNDER the normal CURVE To LEFT OF Z FOR POSITIVE 2. 

Cl 

[1 : 0 5 7 


ARUNCP = PROS (T) /2»C ♦ .500 

oJ 

! 1)072 


DSTaT(I) = (ARUNCP - PCAREA(D) 

OC0077 

7C6 

continue 

c 

“0102 


PRINT 706 

0 

1105 


PRINT 707, (DSTAT(I) ,1=1, NDATA) 

000130 

707 

FORMAT (6(10X,Fl6.5)l 

Of 

7 0130 

708 

FORMAT (//40X,53H 0 VALUES FOR KOLMOGOROv-SMIRNOv GOODNESS OF FIT 




1 TEST/4 IX, 52H (LISTED IN THE SAME ORDER AS CYCLES-TO-FAILURE DATA)/) 

ol 

>3l3c 


return 

000131 


end 


li 

i 

I 

I 

I 

f] 


i 


/ 




t 


C-3. (Cont ’d) . 


cf 0012 

| 

C f 0 012 
0! 0012 

, 1 
1 1 
0 0 

r >oi4 
(' u 0 0 1 5 
C! 7020 

710 

010027 

000030 

r L 

C-- 

0,' 0 037 

c— 

C P 0 0 4 1 
01 1042 

Oy 0044 

711 

000052 

C-- 

0 1053 
000056 
Of ,1062 
0 Jp 0 66 
000102 

712 

olfjl02 

713 

OO0I02 

m 

714 

6 «5 i 02 
066103 



I 

I 

I 

II 


! 


SUBROUTINE ALPHA ( X « NDATA, DATA, XrfF.ANi DEV, ALPHA3* ALPHA4) 
DIMENSION X (NDATA) 

-SUBROUTINE TO CALCULATE THE COEFFICIENTS OF SKt'wNESS AND KURTOsIS 
-CALCULATE THE THIRD MOMENT OF THE DATA (SKEWNESS) 

T0P3 = 0.0 212 

VAR =0.0 

DO 71C I cl, K’DAT A 

VAR = VAR + (X(I) - XMEAN)*»2 

T0P3 = TOP3 *(X(I) - XMEAN>*^3 

SKEW = TOP3 / DATA 

STDEV = SORT (VAR/DATA) 

-ALPHA3 = MOMENT COEFFICIENT OF SKEWNESS. 

ALPHA3 b SKEW/(STDEV««3> 

-CALCULATE THE FOURTH MOMENT OF THE DA.TA (KURTOSIS). 

T0P4 = 0.0 

DO 7U I = 1, NDATA 

TOP4 = TOP4 ♦ C X ( I > - XMEaN)**A 

TKURT = TOP4 / DATA 

-ALPHA4 = MOMENT COEFFICIENT OF KURTOSIS. 

ALPHA4 = TKURT/ (STDEV**4) 

PRINT 712 
PRI N T 713 

PRINT 7 14, ALPHA3, ALPHA4 

FORMAT ( /// 1 9X , 39HM0MENT COEFFICIENT OF SKEV/NESS ( A LPHA3 ) , 1 8X , 39HM 
lOMEf-T COEFFICIENT Of KURTOSIS (ALPHA4)/) 

FORMAT (2l X , 34HF0R NORMAL DISTRIBUTION ALPHA3 a 0,23X,36HFOR NO 
1RMAL DISTRIBUTION ALPHA4 = 3.0/) 

FORMAT (26X,25HFOR ABOVE DATA— ALPHa3 =F6 , 3 , 26X , 25HF0R ABOVE DATA 

1 ALPHA* = , F 6 , 3 ) 

RETURN 

END 


(Cont 'd) . 


- : t~i - i' 
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APPENDIX D 

Program STRENG (FORTRAN) 


D-l. Flow Chart 

D-2. Definition of Variable 


D-3. Computer Listing 


PROGRAM STRENG 





STINT, M,N, RATIO 



D-l. Flow Chart 



215 



NM1 = N-l 

INDEX =1 L = 1 

> 

f 



DSTR = ASTR(J + 1) - ASTR (J) 
DNEAN = AMLOG (J+l) - AMLOG (J) 
DSIG = ASIGL (J + l) - ASIGL (J) 
SLAM (J) = DSTR/ DNEAN 



SLSIG (J) = DSTR/ DSIG 



D-l. Fow Chart (Cont'd). 





216 





217 


0 



f 

PRINT 

(HEADINGS) 


/ 


DOES 
\RATIO = 0 


_NCL_~ 




L. 


YES 


( 

PRINT, STRESS 

{RATIO = INFINITY 


jPRIN T (HEADINGS)^ 




NP - 1 


t 

PRINT ASTR(I), 
ACYCLE(I) , L_ 
AMLOG(I), 
ASIGL(I) 



PRINT 

RATIO 


D-l. Flow Chart (Cont'd) 


















D-l. Flow Chart (Cont'd). 






220 



I NT (4) STR = SMIN 


PRINT 

CHEADINGS) 


DOES " 
\ N = INTZ 


, yNO 

PRINT STR, j 

XMLOG(I) XSIGLC.Il 


PRINT 


(HEADING] 


PRINT J, 
AREA CI.J) 


\ 



STR = STR + SINT(XINT) 4 


D-] . Flow Chart CCont’d). 













"uwwiti «ttnvt*fH(!i( a pw^riwwi iiiwi i wippiiii 















224 



-1. Flow Chart (Cont'd). 



KOLMOGOROV -SMIRNOV GOODNESS OE FIT TEST 


PRINT 

(HEADINGS) 



READ CYTOFR (I,J) 


DOES 

J=NDATA 


DOES > 
T = NP -1 


ANUM(1)=0 
TOTAL = 0 


IF (CYTOFR (I , K)< XCYCLF.(J) 

z 

NDATA = NUM(J) 

THEN AN1 = AN1 + 1 


ANI = 0 







































229 


XJMIN = log cycles of extrapolated minimum value of M 
XJMAX = log cycles of extrapolated maximum value of M 
M = number of cycles (XCYCLES) of XJMIN 

N = number of distributions to be used at a particular stress ratio 

RATIO = stress ratio 

ASTR = bending stress level in psi 

AMLOG = log of bending stress level 

NUM = number of cycles to failure inputs (ie. no. of XCYLCE) 

XCYLCE = the actual number of cycle to failure 

STINT = standard deviation of cycles to failure in log t3rms. 

ACYCLE = number of cycles to failure or mean cycles 

Interpolated sti'ess values 
STR - ASTP 
CYCLE - ACYCLE 
XMLOG - AMLOG 
XSIGL - ASIGL 

AREA (I,J) = cumulative are under log normal cycles to failure curve. 

Parameters of Normal Stress Distribution at Specified Cycles -To-Fai lure . 
CYC = cycles in log value 
X CYCLE - cycles 
SMLOG - mean strength 

SSIGL = the standard deviation of the mean strength 

SIG3M = the minus three sigma limit 

SIG3P = the plus three sigma limit 

SK3 = the skewness value 

SK4 = the value of kurtosis 


D-2. Definitions of Variables. 


ait (; 0 74 
« 00 75 
* 0 07b 
0or.o77 

S OI 

0103 
0001 05 

if 6 107 
*101 11 
*011? 
y.) 6 n s 

#0115 
«1C] 17 

0001?? 

f UC-1 24 
!’ 0 1 20 

00*127 
f^Cl 31 

I >yl32 

t'UO 1 34 
000140 
■ ) 0 1 4 1 
•. .‘Cl 43 
001-144 
00151 
10154 
01*0155 
'00160 
*63 

t',00 I 63 
|j C 1 64 


PROGRAM STRERG ( I NPUT , OUT PUT , T APt ) = I NPUT) 250 

--PROGRAM KINDS NORMAL STRENGTH DISTRIBUTION FROM LOGNORMAL 
•-CTCLES-TG-K AlLURE DISTRIBUTIONS. 

DIMENSION A Si R ( i 0 ) , AMLOGUC) , A3I GL ( 1 0 ) , SLAM (9) , SJG3M (23) » 

1 XMLOG ( 30 1 ) * XSJGL (301 ) • AREA (301 » 23 ) ♦ FREQ (30 1 , 23) , S I G3P (?3 ) * 

2 CYC(23), Sn'LOG (23 ) i SSIGL(23), SK3(23), SK4 ( ?3) ♦ CYTOFR ( 1 0 » 25) , 

3 N U M (10) i INDEX (10) ,OFREG( 10,23) , GFREO ( 10 , 23 ) ,DMAX 0 0,23) »AN?(10) 

4 » CHI SO ( 1 0 « 23 ) .ST (10) »SLS1G(9) ♦ EFREO (10.23) . 1NDF.X2 ( ) 0 ) ,ANUM(lfl> 
5.CYCLES(30l> .ACYCLE12S) . STRLOGUO), STMEAN(23)i XCYClE(23) 

EQUIVALENCE ( AREA , FREO) . 

1 (XS1GL.SSIGL) . (XMLOG, OFREQ) ♦ ( AREA (SCO) t EFREQ) * 

2 (AREniCHlSC) , (AREA (250) iGFREO) 

INI = 1 

i RE AD 1 0 1 x JMI N ♦ X JMA X 1 S 1 NT , M , N . RAT 10 
) FORMAT (2Fl0.6,F10.0,215,F10.5) 

FAILURE DISTRIBUTION VALUES ARE READ IN FROM 
LOWEST TO HIGHEST STRESS AND STRESSES ARE 
•integer values of sint 


1-00024 
(f I 0 0 2 7 


20 

IF (EOF,}) 51 C . 20 

READ 3C. (ASTR(J) ,AHLOG(J) ,ASIGL(J) 

(II 0050 


3C 

FORMAT (FI 0,0 *2F 10.6, 15) 

000050 

*“0063 

3b 


READ 36, (XCYCLE (J) , J=1 ,H) 
FORMAT ( BF 1 0 • 0 ) 

#0063 
6 0 6 0 6 5 

34 


DO 34 1=1, N 

ACYCLE (I)=EXP(AMLOG(l ) ) 


NMl = N-l 

L = 1 • 

INDEX ( 1 ) = 1 

L>0 05 JslfNMl 

05 TR = aSTR (J+ 1 ) -ASTR ( J) 

OMEAN = AMLOG ( J + 1 ) - AfiLOG ( J) 

OSIG = A S I G L ( J + 1 ) - A S I G L ( J ) 

SLAM ( j) = DSTR/D-'-EAN 
IF (OSIG) 40,50.40 
SLSIG(U) = OSTR/ilSIG 
CO r *T I NUt 
JP1 = J*1 
MX = D5TR/SINT 
INDEX (UPl) = L ♦ M X 
SM1N - ASTRO) 

XMLOG(l) = AMLOG(l) 

XSlGL(l) = ASlGLO) 

DO 3 0 Kelt MX 
LX = L ♦ K 

XMLOG1LX) = XMLOG(LX-l ) ♦SlNT/$LAn{ J) 

•1F(US10) 70.60,70 
XSlGL(LX) = XSIGL(LX-l) 

GO TO ftp 

XSIGL(LX) = XS1GL(LX-1)*SINT/SLSIG(J) 

CONTINUE 
L = LX 
. CONTINUE 
PRINT 37 

FORMAT (1H1///44X, 44HNORMAL STRENGTH DISTRIBUTIONS FROM LOGNORMA 
1L/49X, 3?H CYCLES TO FAILURE DISTRIBUTIONS,///) 

IF (RAT 10. EQ. J . 0) GO TO 91 
PRINT d*, RATIO 


D-3. Computer Listing for Program STRENG. 


itt 


if 0172 
1**01 72 

0 0 0 1 7 3 
if; 01 7 7 
(| i ii l 7 7 

01 ' 03 


\h0 
1 bO 


170 

lbo 


FORMAT (44X, 17HEXPERIMENTAL DAT A , BX , UNSTRESS RATIO s,F6.3//l ,,, 

GO TO 93 
PRINT 92 

FORMAT ( 44 X » 1 7 HEX PER 1 MENTAL DAT A , 8X , 23USTRESS RATIO = INFINITY//) 
PRINT 9 4 

FORMAT (BX, 13«STRENGTH HSI . , 9X , j?HMEAN-CYCLES LOG MEAN-CYCLES, 

ldXi 1 1HLOG STD DEV/) 

NP = N-l 

PRINT 93, (ASTR(I), ACYCLF(I), AMLOG ( 1 ) » AS10L ( I ) i Is=2,NP) 

FORMAT {2F20.H,2F20.8) 

PRINT 96 

FORMAT (///54X,26H INTERPOLATED STRESS LEVELS//) 

PRINT 90 

FORMAT (ax,13NSTftEN3TH PS I . , 9X , 32HMEAN-CYCLES LOG MEAN-CYCLES, 

1 8X , 1 1 MLOG STO DE V , I 4x , 1 OHINTEGER (I ) / > 

XI NT > INT 

STR s SHIN 

DOl 10 1*1* LX 1 1 NT 

CYCLES(1)=£xP(XML0G(1) ) ' 

PRINT 1 00 , S TP> , CYCLES ( I ) » XMLOG ( I ) , XSIG.L(I),I 
FORMAT (2F20i?,2F20.D, 120) 

ST K = STR*SINT*XINT 
CONTINUE 
NCLD » N 
N a LX 

CONVERTING LOGNORMAL FAILURE DISI. PARAMETERS AT N 
STRESS LEVELS TO CUMMULATIVE LOGNORMAL FAILURE 
DISTRIBUTION 
DO 1HC’ 1 = 1 »N 
DO iat* J~ 1 , M 
CYC ( J) = ALOG (XCYCLE ( J) ) 

Z~ (CYC (J) -XMLOG (1 ) )/XS16L(I) 

JF (2) 12C ,140,140 
IF (2*4.5) I 60 » 130, 130 
2 a -2 

CALL NORMAL <2,PR08) 

ARt A ( I , J) a ( 1 , 0-PROd ) /2 • 0 
GO TO ISO 

IF (2-3.5) 150,150,170 
CALL NORMAL ( 2 , PR03 ) 

AREA ( I , J) = PROti/2. 0*0.5 
GO. TO l'-*o 
ARE<, ( 1 , J) = 0,0 
GO TO In*} 

ARt A ( 1 , J) = 1 . 0 
CONTINUE 
PRINT 225 
I NT2 - I N T a 4 
STR = SMIN 
PRINT 1*5 

FORMAT (///46X.42HCUMULATI VE LOGNORMAL FAILURE DISTRIBUTIONS) 

00 220 1-1 ,N»1nT2 

PRINT 1 VO, STR, XMLOG (1 ) , XSIGL ( I ) 

FORMAT ( lHf ,5X, liMSTRENGTH = ,F10.0,5X, 9HL0G MEAN , 

1 VriCYCLES = ,F9.6,5X,?OHLOU STD DEVIATION s ,F9.G> 

PRINT 2;o 

FORMAT C 1 X , 3 3RD AT A UELO*' JS J, AND CUMMULATIVE * 

I 13H0IST UP TO J,) 


D-3. (Cont’d) 




232 


\ 


; GfU06 


PRINT 210* (J, AilEAU » J) » JsliH) 

0 uC 1 4 24 

211 

FORMAT (UX,I3,F9.6»)3*F9.6»l3,F9.6.I3iF9.6* 

I 


1 13,F9.6» I3.F9.6, I3.F9.6, I3,F9.b, J3.F9.6) ) 

Oj -424 


STR r. STR4S]NT«XINT*4, 

or, 30 

220 

CONTINUE 

0Uy*33 


00 222 J= 1» MOLD 

Ol l>4 34 


I = INOEX(J) 

04C436 


XNUM = NUM(J) 

0CU437 


DO 224 K = 1 ♦ M 

Of 1-4 41 


OFREQ ( U » K ) = XNUM«4R£A(1 ,K) 

cj ;- ! 4 5 1 

224 

CONTINUE 

, 00U454 

222 

CONTINUE 

! i 

’ C 

FREQUENCY DIST AMD NORMAL DISTRIBUTION PARAMETERS 

0( 1456 


XJM1N = XJM1N-XJ1NT 

o c 0 4 6 0 


PRINT 225 

000464 

225 

FORMAT < 1 Hi ) 

0 1 J4 64 


PRINT 230 

iH J470 

230 

FORMAT (///SIX, 33HCUNULATIVE STRENGTH DISTRIBUTIONS) 

000470 


00 260 J=1 ,M 

Of ) 4 7 2 


XJ = J 

0 1 -.473 


PRINT 240, J* CYC(Ji, XCYCLF. (J) 

GOO DCS 

240 

FORMAT ( 1 HO , 4H J = ,13 , 5X , 1 3HL06 CYCLES = ,F9.6, 

IT 


1 5X « 3hCY CLES = » F 1 0 • 1 ) 

; GID505 


PRINT 250, ( I , AREA ( I , J ) , I = 1 ,N , INT ) 

l 00^523 

250 

FORMAT (?X,3«( I ) , 1X,9HFRE«UENCY/ (IX, I4*F10.6, I5,FI0.6 
1 .15. FI 0. 6, 15, Fie. 6. 15, F10. 6, 15, F1C. 6, 15, FI 0.6, 15, F.l 0.6)) 

[ 0()523 

260 

CON’T ] NUE 

• 01jS26 


NM 1 s N-l 

000530 


DO 27) J= 1 , M 

i?* ; ‘531 


DO 270 I=),NM1 

0| .32 


I PI = 1*1 

000534 


F RtO ( I , J) = AR£A ( I PI , J) -AREA (J , J) 

000542 

270 

CONTINUE 

i 0(1547 


N = NMl 

K 

C 

ME*N AND NORMAL DISTRIBUTION PARAMETERS OF HISTOGRAM 

’ 000 551 


00 30C J= 1 , M 

0JD552 


FI « 0. ' 

^ cIj553 


Flul s't. * 

000554 


U1 = 5M I N~S I NT *0 • 5 

0fejS57 


DO 28) 1=1 ,N 

Op; 560 


UI = UI ♦Sl.NT 

ofo;562 


FI = f 1 *FREO ( 1 , J) 

000566 


FIUI = FIUI*FREQ(I.J)*UI 

0 Pi 5 7 3 

260 

CONTINUE 

010575 


SMUOG(J) = FllJl/Fl 

000577 


.ST ME AN' ( J) = ALOG (S-iLOG ( J) ) 

l 1)603 


SM2 = 0. 

0| .-604 


SM3 = •?. 

OCObC‘5 


SM4 =0.(1 

04*0606 


UI » SM1N“S1NT«0.5 

C -jell 


DO 29C 1 = 1, N 

0 1 0 6 1 2 


UI = U 1 * S J N T 

000614 


SO = (Ul-SMLOG(J) )«(U1~SML0GCJ) ) 

Of 0617 


SM2 = S''.? ♦ SO«FREO{J,J) 

01 624 


5 M3 = S M3 ♦ FHEutl , J)*(UI-SMLO&(J) )«SO 

U(. .32 


SM4 = SM4 4S0*S0« l r REOd ,J) 

Of 0 e 3 7 

290 

CONTINUE 

C 064 2 


Sm2 = SV-2/FI 


I 

1 

i 


D-3. (Cont'd) 


•t 

i 


(r c 0 0 4 3 
( 1 5 0644 
(f 0 6 6 4 5 
000651 
i ! 65 S 
. -60 
0 0 0663 
Q 665 
0670 
000673 
C ''0677 


Ov 0677 
000703 


000703 

<^j 0 7 3 ? 
V 0 6 7 3 5 
0741 
07 4 b 




G 0 0 7 4 5 
(ji )751 

000751 

1 6753 
0755 
r -770 
70 

WO 773 
16774 
0O0775 
[0776 
[0777 

| ooiooi 

I 0^1 002 
1 fill 0 03 
! OV 1 0 1 2 
4 001015 
1 C‘| 1017 

• cl. 1 021 
j 001022 

O ’ 1 0 24 
0 1 0 26 
001032 
(•010 33 
0 1034 
j Oi l 0 37 
1 001041 
11044 



SM3 = 5M3/K1 
$M4 = SM/Fl 
5S1GL(J> = SO’RTF (SM? ) 

SI 03P ( J) = SMLOG(J) ♦ 3» 0*SSIGL ( J) 

$1G3M(J) = SKLOG(J) - 3.0«SSIGL < J) 233 

SK3(J) = SM3/(SSIGt(J)«SM2) 

S*4(J) = SM4/ (SM2«SM2) 

300 CONTINUE 
PRINT 225 
PRINT 305 

305 FORMAT <///3QX|73HPARAMETEN$ OP NORMAL STRESS DISTRIBUTIONS AT SPEC 
1IFIED CYCLES TO FAILURE,//) 

PRINT 310 

• 310 FORMAT ( S 3 X * 20H-3 SIGMA *3 SIGMA /3X » 20HNUM8ER LOG CYCLES , 

1 6X,6MCYCL£S-»10X,13HMEAN STRENGTH, ) OX,9HSTO. DEV. » 7X, 5HL IM IT 1 

2 7 X , 34HL I M I T SKEGNESS KURTOSIS /> 

PRINT 320, ( J,CYC(J) ,XCYCLE(J) iSMLOG(J) » SSIGL(J), 

1 SIG3M ( J) , SI 63P ( J) , SK3 ( J) ,SK4 ( J) , J = 1,M) 

320 FORMAT (4X, I3,F16.6»F13.0»2F20.0»F15.0»Fl2.0i2F12.4) 

PRINT 225 
PRINT 322 
322 FORMAT (1H0) 

C ROUTINE FOR K0LM06UR0 V-SMi «NOV GOODNESS OF FIT TEST 

PRINT 325. 

325 FORMA f (//4 IV i52hO- VALUES FOR KOLMOGOROV-SmJ RnOV GOODNESS OF FIT T 

1EST /) 

DO 600 1=2, NP 
NOATA = NUM(I) 

READ 635, (CYTCFR (1 , J) , J=1 , NDATA) • 

605 FORMAT *(BF10<0) 

600 CONTINUE 

DO 50C- J=1 ,M 
ANUM(1.) e 0,0 
TOTAL = O.C 
DO 665 1=2, NP 
NDATA = NUMU) 

ANl = 0 • 0 

DO 660 K s 1, NDATA 

IF (CYlOFR(l»K) .LE.XCYCLE (J) ) AN1 = AN1 4 l.o 
660 CONTINUE ' ‘ 

AN2(I) = AM 

AN JM ( I ) = ANUM ( I - 1 ) 4 AN1 
TO.TAL = TOTAL 4 ANl 
665 CONTINUE 

00 45!' 1 = 2 » NR 

2 MASTR(J) -SMLOG ( J) )/SSlGL ( J) 

•IF ( Z ) 330,350,350 
330 L - -l 

IF(Z~4*5) 340,36J,360 
340 CALL NORMAL (?,PRO-J) 

PROD A = 0.5-PROB*0.5 
GO TO 37 y 

350 IF (2-4.5) 355065, 3o5 
355 CALL NORMAL (Z,PRCW) 

PROBA - PR0B*0 ,5*0.5 
GO TO 37? 

366 PROD A = 0.0 
GO TO 37 j 


D- 3 . (Cont r d) 


f cl 

: 0 56 

365 

PR03A =1.0 


I poiooo 

370 

6FH£0(li J)=PR0SA 

234 

; c* 

[’064 


IF (TOTAL. £Q .0,3) 60 TO 375 


? i*| 

! 065 


DARLA = ANUMlD/TOTAL 


i 01 

67 


60 TO 380 


! o 

< - -c 6 7 

375 

DARE A = 0.0 


b 

i*0 7 0 

380 

UMAX C I » J) = ASS(DAREA-6FR£Q(I, J) ) 


i 0 

• 1 1 0 1 

46 0 

CONTINUE 


; oul 1 o3 


•PRINT 490, XCYCLE(J), TOTAL , ( AN2 ( I ) 

.OMAX(I.J), I =P, K’P ) 

; C 

125 

490 

FORMAT (6X,F8,0,21H CYCLES TOTAL 

W =,F3.0/(9(F9.0,F6.3) > ) 

i 

125 


PRINT 322 


| 001131 

. 500 

CO r \T I NUc 


i tn 

"134 


60 TO 5 


I O' 

134 

51 C 

STOP 


: 00H36 


END 



I 



l 

l 



D-3. (Cont'd) 



r 

! 


SUuROU T 1 t\E NORMAL (Z,PROt3) 



C 

PP03 “ ThE AREA UNDER NORMAL DISTRIBUTION BETWEEN 



c 

PLUS AND MINUS L STANDARD DEVIATIONS 

G 

! '05 


IF ( 2 - } . 2 ) lOCOi IOOO j 1010 


.i v j 0 7 

1000 

ZSG = 2«Z 

0 

’0010 


PR03 = 0.7978>7455*Z* (0.99999774-ZSG!* (Q, 16659433 


• 


1 -ZSGMO. 024638310 - ZSQ*0 . 0 023974867 )) > 

c 

, 0020 


GO TO 1070 

COCO?} 

1010 

IF (2-2.9) 1020» 1060» 1060 

0 

.'024 

1020 

2 SO = z*z 

0 

• 0 25 


PR08 = 1.0 

000026 


PTERM = 1.0 


'0027 


FACT =1.0 

0 

'027 


000 IN = 3.0 

0 0 0 0 3 1 

1030 

PTERM = -PTERM*ZSQ/(2.0*FACT> 

0 

>0 035 


TERM = PTERM/ODOIN 

0 

j 0 3b 


PROd = PROS + TERM 

0 

.0 037 


IF ( ABS(TERm) - 0.00007} 1 OSO i U)4 0 » 1 04 0 

090043 

104 0 

FACT = FACT ♦ 1.0 

0 

i D 0 4 5 


OODIN = ODDIN +2,0 

C 

[1)047 


GO TO 1 30 

000047 

1050 

PkOB = ).75785455*Z*PR0tj 

°| 

|G 0 5 1 


GO TO 1070 

C 

[JO 52 

10 60 

REC = l./(Z*Z) 

c! 

fO 0 54 


PROB = 1.-0.75755453* EXP (-7XZ/2. 0 ) /Z* 


■to 


1 (1. - PXC*C1. - REC* ( 3 . - REC* ( 1 5 , -REC* 1 05 . ) ) ) ) 

cl 

1)076 

1070 

CONTINUE 

ol 

1)076 : 


RE f URN : 

0 0 ft 0 7 7 


END 


I 

I 

B 

P 
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APPENDIX E 

Listing of Short PDP-8 Programs 


E-l . Program BAR I 
E-2. Program BAR II 
E-3. Program ROTO 

E-4. Program for Least Squares Estimator for Chapter IV Data 
E-5. Program for Least Square Estimator for Chapter V Data 
E-6. Program for Mean Stress Per von Mises-Hencky Ellipse. 
E-7. Program for Slope of Best Fit Equation 


237 


01.10 A "SH" S, "SA" SA; S SM=FSQT(3)*S 


02.10 S C=FSQT(SA*2+SH+2); T %8.1 / "C" C, 
02.20 GT1.1 


* 

* 

* 


"SH" SM, 


i 


v/here : SH = shear stress 

SA = alternating stress ' 

SM “ mean stress 

C = resultant stress vector 
magnitude S^. 

E-l. Listing of POP- 8 Program BAR I. 
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*E A 
*W A 

C-BK HODV 11-219 

*01.10 A "SH" S, "SA" SA; S SM=2*S 

*02.10 S C~PSQT( SA<‘2+SH< , 2 ) ; T %8.1, "C" C # "SK" SH, ! 

*02.20 GT1.1 

*GO 


where j SH = shear stress 

SA - alternating stress 

SM = mean stress 

C - resultant stress vector 
magnitude S r 

E-2. Listing of PDP-8 Program Bar II. 
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* 

*C-8K KODV 11-219 
* 

*01.10 a "R" R, "S" S; S SR=S/FSIH FATN(R) 

iRr 

*02,10 T %8.2 "ROTATED THREE SIGMA" SR, ! 
*02.20 GT 1.10 

**/ 

* 

* 

* 


YJhere R = Stress Ratio 

S - Vertical Strength Distribution’s 
Upper and Lower Three Sigma 
Limits . 


SR - Transformed Strength Distri- 
bution's Upper and Lower Three 
Sigma Limits . 


E-3. Listing of PDP-8 Prgram ROTO. 
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C-FOCAL,1969 

0.10 A "A" A, "3" B, "D" D / "E" E, "H" N 
0.20 S XA=A/E 
01.21 S XB=B/N 

01.23 S YA-D t 2-A + 2 

01.24 S YB=D * 2-B 1* 2 

02.10 S E=XA t 2*YA+XB * 2*YB 

02.11 S J=XA t 4+XBt 4 
02.13 S K=S/J 

0.10 S SU=F3QT(D* 2/K); T "SU" SU! 




Definition of Variables 
A = Sa^L = 

B ~ Sa 2 ~ *2 
D - En 

y a - Sm^ = x^ 
yb = Sta 2 = x 2 



E-4. PDP-8 Program for Least Squares Estimator of the 
Ultimate Strength for Fatigue Data of Chapter IV. 


^C-FOCAL, 1969 

*01.10 a "A" A, "B» B, "C" C , "D" D 
*01.20 S XA=A/3. 5 
*01.21 S XB=B/ .825 
*01.22 S XC~C / .44 
*01.23 S YA=D f 2-A * 2 
*01.24 S YB=D f 2-B 1 2 
*01.25 S YC=D t 2-C t 2 
* 

*02.10 S E=XA* 2*Y'A+XB t 2* YB*XC * 2*YC 
*02.11 S J=XA 4+XBt 4+XCt 4 
*02.13 S K s E/j 
* 

*03.10 S SU=FSQT(D t 2/K); T "SU" SU1 


Definition of Variables 
A “ Sa^ = 

B = Sa 2 = y 2 

C = Sa 3 = y 3 

xa = Sm^ = 

xB = Srug - y ~2 

xC — *• Xg 

E-5 . PDP-8 Program for Least Squares Estimator of 
the Ultimate Strength for Fatigue Data (LSEFD) 
of Chapter V. 


*C-8K MODV 11-219 
* 

*01.10 a "SA" S, ."SB" W; S SH-255300*FSQT(1- <sA/> t Z ) 
*01.20 T 5^.2 "MEAN STRESS" Sii, I 
*01.30 GT 1.10 
** GO 


E-6. 


Listing of PDP-8 Program Which Calculated Mean 
Stress Specified by the Von Hises-lIencJcy Ellip: 



; r 

• ! 



C-8K HODV 11-219 

02.05 A "MO. OF DATA POINTS ",1©,!.;! 

02.10 T “ X - AXIS Y AXIS",! 

02.20 S XS=0 . “ 

02.21 S X0=0 ' 

02.22 S YS=0 

02.23 S YQ=0 

02.24 S XY=0 

02.40 FOR I-l,l,i©; DO 3.0 

02.49 S D=(MD*XQ-XS*XS) 

02.50 S AO=(YS*XQ~XS*XY)/D 
02.60 S Al=(I©*XY-XS*YS)/D 
02.70 S DN=?3QT(D* ( ND*YQ-YS*YS ) ) 

02.75 s r=(md*xy-xs*ys)/d:i 

02.80 T %8.05 ! , "SLOPE ",Al," Y INTERCEPT ",A0, !,! 
02.90 T %6.04 "CORRELATION COEFFICIENT ",R,! 

02.95 Q 

03.10 A " ",X(I)," ",Y(I),I 

03,15 S XS«XS+X(I) 

03.18 S XQ-XO*X(I)*X(I) 

03.20 S YS=YS+Y(I) 

03.24 S YQ- YQ+Y ( I )*Y ( I ) 

03.26 S XY=XY+X ( I ) * Y ( I ) 

*GO 


E-7. Listing of PDP-8 Program which Calculates 
Slope of Best Fit Equation (10.5). 
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